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1. INTRODUCTION 

> . 

Relativistic heavy-ion collisions at center of mass energies in the range of many to 
many hundreds of GeV per nucleon are performed with the ultimate goal of generating 
Q\ \ a new state of matter, the quark-gluon plasma (QGP). To make a QGP, one must heat 
and compress nuclear matter to such high energy densities that the hadrons overlap, 
creating a homogeneous piece of matter consisting only of the stuff usually hidden 
inside hadrons. In the region thus created the quarks and gluons, which are usually 
imprisoned inside the hadrons, become "deconfmed" , i.e. they are able to travel around 
freely over regions which are large compared to the usual confinement length scale of 
about 1 fm (the typical radius of a hadron). Such matter existed in and filled all of the 
volume of the Early Universe during the first microsecond of its life but, as far as we 
know, it has not been recreated in the laborarory anywhere since. 

How do we know our heavy-ion collision experiment has been successful in creating 
QGP? The answer to this question is surprisingly complex. The main reason for this 
is that quarks and gluons cannot travel over large distances outside the hot and dense 
QGP region and thus cannot be detected directly. Everything we can measure in the 
experiment provides therefore only indirect knowledge about the hoped-for QGP state. 
This is in particular true for the bulk of secondary particles produced in the collision 
which are hadrons: they can only be formed towards the end of the collision process 
when the reaction zone, by expansion into the surrounding vacuum, has cooled down 
sufficiently to allow the quarks and gluons to re-hadronize. 

One approach to "prove" the making of QGP (which I stress must still be com- 
plemented by other experimental tests to check for consistency, see Jiirgen Schukraft's 
lectures [ Q for examples) tries to reach a complete understanding of the space-time 
structure and dynamical state of the reaction zone at the "freeze-out point" where the 
measured hadrons decouple. One hopes and by now, from the experiments performed 
at the Brookhaven AGS and CERN SPS during the last decade, has accumulated a 
reasonably convincing body of evidence that at decoupling the matter has reached a 
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state of local thermal and perhaps even chemical equilibrium. It also features strong 
longitudinal and transverse collective (hydrodynamic) expansion. While some of the 
longitudinal expansion may be "primordial", i.e. due to incomplete stopping of the 
two nuclei in the collision region, all of the transverse expansion must have been gen- 
erated dynamically in the reaction. If it is possible to determine the energy density of 
the state at decoupling and simultaneously its expansion velocity, then one can try to 
extrapolate this state backwards in time to a point of vanishing transverse expansion 
and check whether there the energy density was above the critical value of about 1 
GeV/fm 3 where one expects the phase transition to the QGP to occur. 

But how does one measure the energy density at decoupling? To measure the 
total energy of the state is not too difficult: one measures the momenta of all emitted 
particles and their masses and adds them up. But to compute the energy density we 
must also know the volume of the reaction zone. And there is no known way to measure 
this volume directly. The collision fireball just doesn't hang around long enough so that 
we could shine light on it (or a suitable other species of particles) and measure its size 
by diffraction. 

The only known way to obtain indirect experimental information on the space- 
time structure of the particle emitting source created in a relativistic nuclear collision is 
through two-particle intensity (Hanbury-Brown-Twiss (HBT)) interferometry [[J. The 
goal of this method is to extract the space-time structure of the source from momentum 
spectra which are the only measurable quantities, making use of the quantum statistical 
correlations between pairs of identical particles. The basic idea, presented in a very 
naive and oversimplified way which to correct I will spend three hours of lectures 
during this Advanced Summer Institute, is as follows: Consider a source (for simplicity 
spherically symmetric) with radius R, and the emission of a pair of identical particles 
from point X\ in the source with momentum pi and point x 2 with momentum p 2 . If 
these two points are well separated in phase-space, i.e. they satisfy 

(x\ - 4)(Pi _ Pi) > > i = lori = 2ori = 3, (1) 

this process can be treated classically. If on the other hand, 

{x[ - xi)(p\ - pi) <27rh, i = 1, 2 and 3 , (2) 

the two particles sit close in phase-space, and quantum mechanics can no longer be 
ignored. The most important quantum mechanical correction to be taken into account 
is the (anti-)symmetrization of the two-particle wave function: it ensures vanishing 
probability for two identical fermions to originate from the same phase-space point, 
and for bosons it leads to an enhanced probability to find them at the same point in 
phase-space compared to the classical expectation (bosons are "social subjects"). 

Since the distance in coordinate space (x\ — x l 2 ) is limited by the finite diameter 2R 
of the source, we can force the system into the quantum domain by measuring particle 
pairs with smaller and smaller relative momentum p\—pi- Once q% = p\ —p\ becomes 
smaller than irh/R, the two particles can no longer avoid quantum mechanics by escap- 
ing to larger relative distances in coordinate space. There emission probability will be 
affected by wave function symmetrization, leading in the case of bosons (fermions) to 
an enhanced (reduced) pair emission probability compared to the classical expectation 
(which would be simply the product of the individual single-particle emission proba- 
bilities). The two-particle correlation function is thus expected to begin to appreciably 
deviate from unity for relative momenta q < q* ~ (%/R). The critical value q* at which 
this effect sets in (conveniently one chooses for q* the value where the correlation is half 



way between it maximum (minimum) value and 1) is thus a measure for the geometric 
radius R of the source. 

So far the naive picture of how two-particle HBT interferometry works. Unfor- 
tunately, the only situation where it applies more or less directly is for photon in- 
terferometry of stars for which the method was invented. I will spend my time here 
in Dronten to explain to you why this naive picture is generically wrong and how to 
substitute it correctly. The basic reason why the above simple cartoon ceases to work 
in high energy nuclear and particle physics is that the sources created in hadronic or 
heavy-ion collisions live only for very short time periods and feature inhomogeneous 
temperature profiles and strong collective dynamical expansion. I will show you that 
for such sources the HBT radius parameters (half widths of the correlation function) 
generally don't measure the full source size, but only so-called "space-time regions of 
homogeneity" inside which the momentum distribution varies sufficiently little so that 
the particles can actually show the quantum statistical correlations. They also mix the 
spatial and temporal structure of the source, and we will learn tricks (in particular a 
new way of parametrizing the correlation functions) which are designed unfold these 
different aspects from the data. 

The size of the just mentioned homogeneity regions varies with the momentum 
of the emitted particles, causing an important dependence of the HBT parameters on 
the pair momentum. I will show you how this momentum dependence can be used to 
extract the strength of the collective flow of the source at decoupling. To do so in a 
quantitative way requires the detailed consideration of many physical features of the 
particle emission process. One such feature which I will discuss here in some detail is 
the fact that often some of the particles one uses in constructing the pair-correlation 
function don't come directly from the source but are created well after decoupling by the 
decay of unstable resonances. Resonance decays not only affect the size and momentum 
dependence of some of the HBT radius parameters but also the overall strength of the 
correlation as described by the so-called "chaoticity parameter" A. 

My lectures are structured in the following way: In the First Lecture, I discuss 
the general connection of the measured one- and two-particle spectra with the phase- 
space distribution of particles in the source. In particular I discuss the important 
aspects of "chaotic" versus "coherent" particle emission and how to implement them 
in a formal way. This establishes the formal background which I exploit in the Second 
Lecture to extract in a more quantitative way general ( "model-independent" ) analytic 
relationships between the geometric and dynamic space-time structure of the source 
and certain features of the pair correlation function. In the Third Lecture I analyze 
these relationships quantitatively in the framework of a general class of model sources 
with finite geometric extension and three-dimensional collective expansion. The results 
should give you a good qualitative and even semi-quantitative feeling of the expected 
behaviour of the pair correlation function in relativistic heavy-ion collisions. 

The Fourth Lecture should have presented a comparison of the calculations with the 
data, but fortunately the organizers gave me only three hours time so that I was spared 
the embarrassment of having to admit that at the present moment no quantitative such 
comparison exists: the theoretical analysis is still so new that most of the diagrams 
you will see are less than a few months old, and the new high-quality heavy-ion data 
which allow for the detailed multidimensional analysis advocated here are still so hot 
("preliminary") that the experimentalists wouldn't give me permission to show them 
anyhow. Please look into Barbara Jacak's lecture notes in this volume [ [IJ and into the 
Proceedings of the Quark Matter '96 Conference in Heidelberg, 20-24 May 1996, if you 
want to catch a first glimpse of the data. 



2. LECTURE 1: SPECTRA AND EMISSION FUNCTION 



2.1. One- and two-particle spectra 

The covariant single- and two-particle distributions are denned by 

dN 

Pxip) = E— = E(a;a p ), (3) 

P2(p a ,Pb) = E a E b ' N = E a E b (a+ a+ a Pb a Pa ) , (4) 

where a+ (a p ) creates (destroys) a particle with momentum p. The angular brackets 
denote an ensemble average, 

(6) = tr (p6) , (5) 

where p is the density operator associated with the ensemble. (When talking about an 
ensemble we may think of either a single large, thermalized source, or a large number of 
similar, but not necessarily thermalized collision events.) The single-particle spectrum 
is normalized to the average number of particles, (N), per collision, 

^PM = (N), (6) 

while the two-particle distribution is normalized to the number of possible pairs, 
(iV(iV - 1)), per event: 



d 3 pa (Ppp 

The two-particle correlation function is defined as 



f^ d ^P2( Pa ,P b ) = (N(N-l)). (7) 



C{Pa ' Pb) - (N(N-l)) PMPip,)- (8) 

If the two particles are emitted independently and final state interactions are neglected 
I will show that it is possible to prove a generalized Wick theorem, 

P2(Pa, Pb ) = {N( ^ N ~ (Pl(p a )Px(p h ) ± \S(p a , Pb )\ 2 ) , (9) 

where the plus (minus) sign refers to bosons (fermions), and we have defined the fol- 
lowing covariant quantity: 



S(p a ,Pb) = \jE a E b (a+ a Pb ) . (10) 

If the ensemble corresponds to a thermalized source in global thermodynamic equilib- 
rium, i.e. p is the (grand) canonical density operator, this is just the well-known "ther- 
mal Wick theorem" [ || . But even in nonthermal sources emission may be "chaotic" , 
i.e. the emission of particle a at point x a may be completely independent from the 
emission of particle b at x b . This can be described by assigning the wave function of 
the emitted particle a random phase <fi which is averaged over in the ensemble average. 
I will show that this is also sufficient to guarantee Wick's theorem (|10D . - The opposite 
to independent, "chaotic" particle emission is emission from a coherent state (e.g. in a 
laser) where the phases of all emitted particles are fully correlated with each other; in 



that case the second term in (|T0| ) is completely missing. The general situation can be 
parametrized by a superposition of density operators, 



P = « Pchaotic + (1 - ®) Pcoherent , (H) 

with < a < 1, where only the chaotic first part contributes to the second term in (|9|): 

(a+a£a b a a ) = (a+a a ) (a£a b ) ± a\ (a+a b ) ch \ 2 

+ a(l - a) ((a+a a ) ch - (a^a a )coh) ((ata b ) ch - (a^a&) cc .h) . (12) 

The last term contributes only if the chaotic and coherent parts of the density operator 
generate different single-particle spectra. One easily checks that for a < 1 in such an 
ensemble the correlations are "incomplete" , i.e. at p a = p b the correlation function (§) 
approaches a value less than 2 for bosons and larger than for fermions. - Since so 
far the heavy ion data don't appear to require a coherent component in the particle 
production process, I will for the remainder of these lectures assume a = 1 and refer 
the interested reader for the possible effects from partial coherence to the literature [ 

Assuming the generalized Wick theorem (^) the correlation function (|8|) can be 
written as 

( a p a a Pa}(% b a Pb ) 

Note that the second term is positive definite, i.e. the correlation function cannot, for 
example, oscillate around unity. [If you see such a behaviour in the literature [ § (and 
the authors did not include final state interactions) it is wrong.] 

From now on I will assume that the emitted particles are bosons, and for con- 
venience I will call them pions, although nearly everything in the first two lectures 
applies equally well to other bosonic particles (except where explicitly stated). In the 
last Lecture I will be a little more specific and distinguish between 2-pion and 2-kaon 
correlations for comparison. 



2.2. The generalized Wick theorem 



Pions are created in heavy-ion collisions throughout the history of the collision 
process, but we are only interested in those pions which reach the detector as free, 
non-interacting particles. (Unfortunately, pion interferometry can only be practically 
performed with charged pions, because neutral pions decay too rapidly. But charged 
pions have a long-range final state interaction, the Coulomb repulsion from the other 
pion in the pair with equal charge, as well as the attractive or repulsive Coulomb 
interaction with the hundreds of other produced charged particles, including the charge 
of the protons in the fireball. A proper treatment of this many-body Coulomb problem 
is a difficult task [ |J . Since I don't know yet how to do it properly, I will simply neglect 
the Coulomb interaction in the final state - assuming that somehow the experimentalists 
know how to approximately correct their measured correlations for it [ |3|.) 

We assume that the last interaction of the pion, in which the finally observed 
pion is created in its free asymptotic state, can be parametrized by a classical source 
amplitude J(x). The solution of the free Klein-Gordon equation for pions generated by 
such a classical current, 

(□ + m 2 )0(x) = J(x) , (14) 
with outgoing boundary conditions is given [ || by a classical ( "coherent" ) state 

| J) = e""/ 2 exp (i J d 3 P J(p) a+) |0> (15) 



where 

J(p) = / . X exp[i(E p t - p-x)] J{x) (16) 



\2-Kf2Ep 

is the on-shell Fourier transform of the source J(x), and the normalization of the state 
is given by 

n= fd 3 p\J(p)\ 2 . (17) 



The state flT5|) is an eigenstate of the destruction operator: 

a p \J)=iJ{p)\J). (18) 

2.2.1. Emission from a single coherent state 

If there is only a single classical source J(x), the corresponding density operator of 
the "ensemble" is just the projection operator on the coherent state, ^coherent = \J)(J\, 
and, using ([18|) , the single- and two-particle spectra (Hd) are easily evaluated: 

E^- = E{J\a1 p a p \J) = E\J(p)\\ (19) 

EaEb d*p N d 3 Pb = E * E b(JK«K b <<\ J ) = E *\J(p*)\ 2 - E b\J(Pb)\ 2 - (20) 

Obviously, there is no exchange term in the two-particle spectrum which is simply given 
by the product of the single-particle spectra. A coherent state thus has no Bose-Einstein 
correlations. 

2.2.2. Emission by a chaotic superposition of classical sources 

This changes if we consider a superposition of classical source amplitudes each of 
which emits free pions independently, i.e. with a random phase [ ||: 

N 

J(x) = £ e ** e~ ipi < x - x J J (x - Xi) . (21) 
i=i 



The construction rule [ |10|, 11] for this source is obvious: we take N sources Jq(x) 
with identical internal structure, give each of them a boost with 4-momentum then 
translate them to different positions Xi in the fireball and supply them with a random 
phase (pi. This allows for arbitrary x-p correlations [ |10[ (i.e. correlations between 
the momentum spectrum of the emitted particles and the point from where they are 
emitted). The momenta Pi of the sources can, but need not be on the pion mass-shell; 
for example, the source could be a decaying A-resonance with 3-momentum pi. The 
on-shell Fourier transform of fl2"T| ) is 

N 

J(p) = Y,e^e^>J (p- Pl ), (22) 
i=i 

where 

Up -i*)= f r, d =k Sp ~ p&x j ^ ( 23 ) 

'{2n)*2E p 



is the (regular) Fourier transform of Jq(x), and p is on-shell while Pi may be off-shell. 
The state | J) which is defined by inserting (p^) into ( |15|) now depends on the parameters 
{xi,pi,4>i;i = 1,...,N}: 

\J) = \J[N;{x,p,<f>})). (24) 

The ensemble of sources can be defined in terms of a density operator p which fixes 
the distribution of these parameters. We assume that the number of sources N is 
distributed with a probability distribution P N , the phases are distributed randomly 
between and 2n, and the source positions Xi and momenta Pi are distributed with a 
phase-space density p(x,p), with normalizations 

oo oo „ 

"£P N = 1, Y, NP N = (N), d 4 xd 4 pp(x,p) = l. (25) 

N=0 N=0 J 

The corresponding ensemble average is given by 

tr(pO)= E^Il d 4 x t d 4 p t p(x t , Pl ) / -^{J[N',{x,p,<i>}]\0\J[N-,{x,p,<j>}]). 
n=o l= i J Jo 2n 

(26) 

The calculation of the single-particle spectrum is straightforward: 

00 N r r 2n d(f>i 

(a+a p ) = ^PnYI d 4 x i d 4 p i p(x i ,p i ) - 1 

n=o i=l J Jo Z7X 

N 

x £ e^-^e i ^ x "- x ^J*(p-p n ,)J (p-p n ). (27) 

n,n'=l 

After performing the integrations over the phases <pi in the double sum over n and n', 
only the diagonal terms with n=n' survive. For each term in the remaining single sum 
over n the integrations over Xj and pi, i ^ n, can be done using the normalization 
condition (|25|) . After suitably relabelling the dummy integration variables for the one 
remaining x- and p-integration we end up with N identical terms under the sum over 
n. This allows to perform the sum over N, and we simply get 

P 1 (p) = E p (\J(p)\ 2 ) = (N)E p fd 4 x'd 4 p'p(x',p')\J (p-p')\ 2 (28) 

= (N) E p [ d 4 p> p(p') \J (p - p>)\ 2 . (29) 



The single particle spectrum is thus obtained by folding the intrinsic momentum spec- 
trum | Jo (p) | 2 of the individual source currents Jo with the 4-momentum distribution of 
the sources, p(p) = J d 4 xp(x,p). 

The algebra for the two-particle spectrum is a little more involved. It is useful to 
first compute 

00 N f f 27T d(pi 

(Ka^pt) = J2PnJ[ d 4 x i d 4 p i p( y x i ,p i ) 

7V=0 i=\ J Jo Z7T 

N 

£ jitn-K') e i(W*n-p..* n ,) J*(p a _ p n ,)J (p b - p n ) . (30) 



N 

X 

n,n'=l 



Again only the terms n—n' survive the phase average, and after doing the dummy 
integrations over Xi,pi, i ^ n, one finds that the remaining sum over n contains again 
iV identical terms, such that the sum over N can be performed: 

(K a O = ( N ) [ d 4 x'd 4 p'p(x',p')e^'^ x ' J*( Pa -p')J (p b -p'). (31) 



With this auxiliary result at hand we can now attack the two-particle spectrum. From 
the definitions one finds 



oo N . „ 2 7T 



N=0 i=l ' 
N 

x E 

n,n'm,m'=l 



o 2tt 



X Jq (Pa ~ Pn')Jo(Pb ~ Pm')Jo(Pb ~ Pm)Jo(Pa ~ Pn) ■ (32) 

The integration over the phases 0, now yields two types of nonvanishing contributions: 
n = n',m = m! and n = m', m = n'. The term where all four summation indices are 
equal, n = m = n' = m', should be omitted [ [5|: it corresponds to emission of both 
particles from the same elementary source, and if one carefully first puts the whole 
system in a finite volume V, performs the calculation there and lets V — > oo in the 
end, then this term is suppressed relative to the others by a factor 1/V. We thus get 



N 



N 



= E -^11 / d i Xid i p i p{Xi,p i ) E [iMPa ~Pn)\ 2 \Jo(Pb ~ Pm)\' 
N=0 i=l n^m 

+ e l{Pa ~ PbHXn - Xm) J*( Pa - Pm)Jo(Pb - Pn)J (Pb - Pm)Jo(Pa ~ Pn) 
N 



1 Y n 

E ^11 / d^x i d' i p i p{xi,p i 



N=0 i=l 
N 



N 



E IMPa ~Pn)\ 2 E \ J o(Pb-Pr, 



N 



\m=l 



+ E J^ Pb -Pn)J (Pa-Pn) 

X (E e- i( **-P»>*™ J (p b -Pm)JZ(pa-Pm)) ■ 
\m=l / . 



(33) 



After again doing the dummy integrations over Xi,pi, i ^ n,m, one realizes that each 
of the two terms in the square bracket contains N(N — 1) identical terms, yielding a 
factor (N(N — 1)) after performing the sum over N. Up to this factor, the first term 
is just a product of two terms of the type (|28f) , i.e. a product of single-particle spectra, 
while the second term is recognized as a product of fl3~T|) and its complex conjugate. 
We thus have 



P2(Pa,Pb) 



(N(N-l)) 

(AO 2 
(N(N-l)) 

(AO 2 



E a E b 
E a E b 



(«p a « P J(«p 6 + \K a & Pb)t 



(34) 



(\APa)\ 2 )(\APb)\ 2 ) + \(r(Pa)J(Pb))\ 2 , (35) 



which proves the generalized Wick theorem (^). 

2.3. Source Wigner function and spectra 

These expressions can be rewritten in a very nice and suggestive way by introducing 
the so-called "emission function" S(x, K) [ [12], [13], 0: 



S(x,K) 



d A y 
2(2tt) ; 



(36) 



It is the Wigner transform of the density matrix associated with the classical source 
amplitudes J(x). This Wigner density is a quantum mechanical object defined in phase- 
space (x, K); in general it is neither positive definite nor real. But, when integrated over 
x or K it yields the classical (positive definite and real) source density in momentum or 
coordinate space, respectively, in exactly the same way as a classical phase-space density 
would behave. Furthermore, textbooks on Wigner functions show that their non-reality 
and non-positivity are genuine quantum effects resulting from the uncertainty relation 
and are concentrated at short phase-space distances; when the Wigner function is 
averaged over phase-space volumes which are large compared to the volume (27ih) 3 
of an elementary phase-space cell, the result is real and positive definite and behaves 
exactly like a classical phase-space density. 

The emission function S(x, K) is thus the quantum mechanical analogue of the 
classical phase-space distribution which gives the probability of finding at point x a 
source which emits free pions with momentum K. Please note that K in S(x, K) can 
be off-shell. Also, it is defined in terms of a 4-dimensional Wigner transform of the 
source density matrix [ O, in contrast to the 3-dimensional expression suggested by 



Pratt [ [T3|] which neglects retardation and off-shell effects. 

Using Eq. flTB|) it is easy to establish the following relationship: 

d A xi d 4 X2 



J*(Pa)J(Pb) = J j—^j==exp(-ip a -x 1 + ip b -x 2 )J*{x 1 )J(x 2 ) 

exp(—iq-x — iK-y)J*(x + \y)J(x — \y) , (37) 



(2n) 3 2^KE b 
d 4 x d A y 

where x = |(xix 2 ) and y = x\ — x 2 . Inserting this into Eqs. ( ^8f) and (|35|) one finds 
the fundamental relations: 

E K §^ = Jd 4 xS(x,K), (38) 

C(aK) = 1 | \Sd*xS(x,K)e**\ 2 
{q ' ' ~ + Jd^xS(x,K + ^q) Jd 4 xS(x,K-y)- { } 



For the single-particle spectrum (|38|) , the Wigner function S(x,K) on the r.h.s. must 
be evaluated on-shell, i.e. at K° = Ek = V m 2 + K 2 . For the correlator (|39|) we have 
defined the relative momentum q = p a — pb, q° = E a — Ef> between the two particles 
in the pair, and the total momentum of the pair K = (p a + p&)/2, K° = (E a + E h )/2. 
Of course, since the 4-momenta p a & of the two measured particles are on-shell, p° = 



Ei = ym 2 + p 2 , the 4-momenta q and K are in general off-shell. They satisfy the 
orthogonality relation 

q-K = 0. (40) 

Thus, the Wigner function on the r.h.s. of Eq. (|39D is not evaluated at the on-shell 
point K° = Ek- This implies that for the correlator, in principle, we need to know the 
off-shell behaviour of the emission function, i.e. the quantum mechanical structure of 
the source. Obviously, this makes the problem appear rather untractable! 

Fortunately, nature is nice to us: the interesting behaviour of the correlator (its 
deviation from unity) is concentrated at small values of \q\. Expanding K° = (E a + 
E\,)/ 2 for small q one finds 

/f °^( 1 + 4 + (t)H A " (41) 



Since the relevant range of q is given by the inverse size of the source (more properly: the 
inverse size of the regions of homogeneity in the source - see Lecture 2), the validity 
of this approximation is ensured in practice as long as the Compton wavelength of 
the particles is small compared to this "source size". For the case of pion, kaon, or 
proton interferometry for heavy-ion collisions this is true automatically due to the 
rest mass of the particles: even for pions at rest, the Compton wavelength of 1.4 
fm is comfortably smaller than any typical nuclear source size. This is of enormous 
practical importance because it allows you essentially to replace the source Wigner 
density by a classical phase-space distribution function for on-shell particles. This 
provides a necessary theoretical foundation for the calculation of HBT correlations from 
classical hydrodynamic or kinetic (e.g. cascade or molecular dynamics) simulations of 
the collision. 

In photon interferometry there is no rest mass available to help you: for photons, 
the approximation K° « Ek can only be justified if they escape from the source with 
high momentum, and in HBT interferometry with soft photons the quantum mechanical 
nature of the emission function needs to be explicitly considered. In practice this means 
that one must study the photon production processes microscopically and quantum 
mechanically. 

If the single-particle spectrum is an exponential function of the energy then it is 
easy to prove [ 14] that one can replace the product of single-particle distributions in 



the denominator of (p9|) by the square of the single-particle spectrum evaluated at the 
average momentum K: 



J d A x e iq -' x S(x, K) 



Jd 4 xS(x,K) 



1 



2 



(e lQ - x ) . (42) 



The deviations from this approximation are proportional to the curvature of the single- 
particle distribution in logarithmic representation [ [14| . They are small in practice 
because the measured single-particle spectra are usually more or less exponential. In 
the second equality of fl4"2] ) we defined (. . .) as the average taken with the emission 
function; due to the independence of S(x, K) this average is a function of K. This 
notation will be used extensively in Lecture 2. 

The ensemble average on the r.h.s. of (|36|) is defined in the sense of Eq. (^6j) and 
can be evaluated with the help of the definition fl2"T|). One finds 

S(x, K) = (N) J d 4 z d 4 q p(x - z, q) S (z, K - q) , (43) 

where 

So(x, P ) = J J^y e -**J5(x + ly)J {x-\y) (44) 

is the Wigner function associated with an individual source Jo- This establishes a 
similar folding relation for the Wigner function itself as we have already obtained in 
(p9f) for the single-particle spectrum: the emission function of the complete source 
is obtained by folding the Wigner function for an individual pion source Jo with the 
Wigner distribution p of these sources. Eq. ( |43"D is useful for the calculation of quan- 
tum statistical correlations from classical Monte Carlo event generators for heavy-ion 
collisions: (N)p(x,p) can be considered as the distribution of the classical phase-space 
coordinates of the pion emitters (decaying resonances or 2-body collision systems), and 
S (x,p) as the Wigner function of the free pions emitted at these points (for example, 
a Gaussian in Quantum Molecular Dynamics calculations [ [TB|). Replacing the former 
by a sum of ^-functions describing the space-time locations of the last interactions and 



the pion momenta just afterwards, and the latter by a product of two Gaussians with 
momentum spread Ap and coordinate spread Arc such that AxAp > h/2, we recover 



the expressions derived in [ 16 



The fundamental relations fl38|) and (^) resp. ( f42|) show that both the single- 
particle spectrum and the two-particle correlation function can be expressed as simple 
integrals over the emission function. The emission function thus is the crucial ingredient 
in the theory of HBT interferometry: if it is known, the calculation of one- and two- 
particle spectra is straightforward (even if the evaluation of the integrals may in some 
cases be technically involved); more interestingly, measurements of the one- and two- 
particle spectra provide access to the emission function and thus to the space-time 
structure of the source. This latter aspect is, of course, the motivation for exploiting 
HBT in practice. In my second and third Lecture I will concentrate on the question to 
what extent this access to the space-time structure from only momentum-space data 
really works, whether it is complete, and (since we will find it is not and HBT analyses 
will thus be necessarily model-dependent) what can be reliably said about the extension 
and dynamical space-time structure of the source anyhow, based on a minimal set of 
intuitive and highly suggestive model assumptions. 



3. LECTURE 2: MODEL-INDEPENDENT DISCUSSION OF HBT 
CORRELATION FUNCTIONS 

In this lecture I will discuss very general relations between the space-time structure 
of the source (as encoded in the x-dependence of the emission function S(x, K)) and the 
shape of the two-particle correlation function. These relations are valid for arbitrary 
emission functions, and in this sense the discussion is model-independent. It never- 
theless provides important insight into the physical features of HBT interferometry, in 
particular for short-lived dynamical sources, and it clarifies what HBT can achieve and 
what not. 

3.1. The mass-shell constraint 

Expressions (^,^) show that the correlation function is related to the emission 
function by a Fourier transformation. At first sight this might suggest that one should 
easily be able to reconstruct the emission function from the measured correlation func- 
tion by inverse Fourier transformation, the single particle spectrum (^) providing the 
normalization. This is, however, not correct. The reason is that, since the correlation 
function is constructed from the on-shell momenta of the measured particle pairs, not 
all four components of the relative momentum q occurring on the r.h.s. of (^) are inde- 
pendent. They are related by the "mass-shell constraint" (|4*0|) which can, for instance, 
be solved for q°: 

q° = f3-q with (3 = ^^^-. (45) 

(3 is (approximately) the velocity of the cm. of the particle pair. The Fourier transform 
in ( [4"2p is therefore not invertible, and the reconstruction of the space-time structure of 
the source from HBT measurements will thus always require additional model assump- 
tions. 



It is instructive to insert ( |45| ) into (f42|): 

/ d 4 x exp(iq-(x — fit)) S(x, K) 



C(q,K)^l + 



fd A xS(x,K) 



(46) 



This shows that the correlator C(q,K) actually mixes the spatial and temporal in- 
formation on the source in a non-trivial way which depends on the pair velocity f3. 
Only for time-independent sources things seem to be simple: the correlator then just 
measures the Fourier transform of the spatial source distribution. Closer inspection 
shows, however, that it does so only in the directions perpendicular to (3 since the time 
integration leads to a (^-function 6~((3-q): 



lim 



/5t dt exp (— i qf3 1) 



I-T dt 



hm ^ 6(q-(3) . (47) 



This implies that there are no correlations in the direction parallel to the pair velocity (3 
(which will be called the "outward" direction below), i.e. C = 1 for g out ^ 0. The width 
of the correlator in this direction vanishes! This should puzzle you: wouldn't you have 
thought that the width of the correlator in the "outward" direction is inversely related 
to the source size in that direction (which is, of course, perfectly finite)? As we will 
see in the next subsection this unexpected behaviour is just another consequence of the 
mixing of the spatial and temporal structure of the source in the correlator: The width 
parameter of the correlator in the "outward" direction receives also a contribution from 
the lifetime of the source which in this case diverges, leading to the vanishing width of 
the correlator. 



3.2. ^-dependence of the correlator 

Eq. (f42|) shows that in general the correlator is a function of both q and K. Only if 
the emission function factorizes in x and K, S(x, K) = F(x) G(K), which means that 
every point x in the source emits particles with the same momentum spectrum G(K) 
(no "x-fT-correlations"), the .fT- dependence in G(K) cancels between numerator and 
denominator of (0), and the correlator seems to be .fT-independent. However, not even 
this is really true: even after the cancellation of the explicit K- dependence G(K), there 
remains an implicit independence via the pair velocity (3 ~ K / Ek in the exponent on 
the r.h.s. of Eq. (f46|) ! Only if both conditions, factorization of the emission function 
in x and K and time-independence of the source, apply simultaneously, the correlation 
function is truely i^-independent (because then the /3-dependence resides only in the 
5-function fl47|)). 

The only practical situation which I know where this occurs and a .ff-independent 
correlation function should thus be expected is in HBT interferometry of stars for 
which the method was invented [ |T7[. It is hard to believe that this complication in the 
application of the original HBT idea to high-energy collisions went nearly unnoticed for 
more than 20 years and was stumbled upon more or less empirically by Scott Pratt in 
his pioneering work on HBT interferometry for heavy- ion collisions [ |13j only in 1984! 

If one parametrises it by a Gaussian in q (see below) this means that in general 
the parameters ( "HBT radii" ) depend on K . Typical sources of x-K correlations in the 
emission function are a collective expansion of the emitter and/or temperature gradients 
in the particle source: in both cases the momentum spectrum ~ exp[— p-u(x)/T(x)} of 
the emitted particles (where u^(x) is the 4- velocity of the expansion flow) depends 
on the emission point. In the case of collective expansion, the spectra from different 
emission points are Doppler shifted relative to each other. If there are temperature 
gradients, e.g. a high temperature in the center and cooler matter at the edges, the 
source will look smaller for high-momentum particles (which come mostly from the hot 
center) than for low-momentum ones (which receive larger contributions also from the 
cooler outward regions). 



We thus see that collective expansion of the source induces a independence of the 
correlation function. But so do temperature gradients. The crucial question is: does 
a careful measurement of the correlation function, in particular of its independence, 
permit a separation of such effects, i.e. can the collective dynamics of the source be 
quantitatively determined through HBT experiments? We will see that this is not an 
easy task; however, with sufficiently good data, it should be possible. In any case, the 
independence of the correlator is a decisive feature which puts the HBT game into a 
completely new ball park. Even if it sounds exaggerated and may at first offend some 
of my experimentalist friends who are busy fighting the limited statistics of their data: 
two-particle correlation measurements which are not able to resolve the independence 
of the HBT parameters are, in high energy nuclear and particle physics, essentially 
useless. [Unfortunately, this applies to all the HBT data from pp and e + e~ collisions 
which I am aware of. In my opinion, a renewed investigation of two-particle correlations 
from pp and e + e~ collisions, using the powerful new tool of multidimensional HBT 
analysis, should be a high priority project - as it is, we have practically nothing with 
which to compare our heavy-ion results in a meaningful way] 

3.3. The Gaussian approximation 



As motivated in the Introduction, the most interesting feature of the two-particle 
correlation function is its half-width. Actually, since the relative momentum q = P1—P2 
has three Cartesian components, the fall-off of the correlator for increasing q is not 



described by a single half-width, but rather by a (symmetric) 3x3 tensor [ 18 1 which 
describes the curvature of the correlation function near q = 0. We will see that in 
fact nearly all relevant information that can be extracted from the correlation function 
resides in the 6 independent components of this tensor [ 19]. This in turn implies 
that in order to compute the correlation function C it is sufficient to approximate the 
source function S by a Gaussian in x which contains only information on its space-time 
moments up to second order. [Gaussian approximations for the emission function have 
been used for the discussion of HBT correlation functions in many different variants [ 
14 , lBl |19| , p0| , 2lj |22|, [23, 24f, a perfect example how research proceeds by trial and 
error. Here I give the rigorous derivation first published in [ Q .] 



Let us write the arbitrary emission function S(x, K) in the following form: 



S(x,K) = N(K) S{x{K),K) e 



\x^{K)B^{K)x v {K) 



+ 6S(x, K) , (48) 



where we adjust the parameters N(K), x^(K), and B^ V (K) of the Gaussian first term 
in such a way that the correction term SS has vanishing zeroth, first and second order 
space-time moments: 



J d 4 x5S{x,K) = J d*xx»5S{x,K) = J d 4 x x^x v 5S(x, K) = 0. 



This is achieved by choosing 

N(K) 

(b- 1 ) (K) 

\ i fit/ 



E 



dN det B^K) 



K " 



d 3 K S(x(K),K) 



(49) 



(50) 

(51) 
(52) 



The (independent) average over the source function (. . .) has been defined in Eq. (|412|). 
The normalization factor ( |50|) ensures that the Gaussian term in ( 48 ) gives the correct 



single-particle spectrum ([38]); it fixes the normalization on-shell, i.e. for K° = Ex, but 
as we discussed this is where we need the emission function also for the computation 



of the correlator. (Note that for photon interferometry this may not be true, and ([50]) 



should then be replaced by a suitable generalization.) x(K) in ( 51) is the centre of the 



emission function S(x, K) and approximately equal to its "saddle point", i.e. the point 
of highest emissivity for particles with momentum K. The second equality in (|52|) de- 
fines x as the space-time coordinate relative to the centre of the emission function; only 
this quantity enters the further discussion, since, due to the invariance of the momen- 
tum spectra under arbitrary translations of the source in coordinate space, the absolute 
position of the emission point is not measurable in experiments which determine only 
particle momenta. Since x(K) is not measurable, neither is the normalization N(K) [ 
p3| as its definition (|50| ) involves the emission function at x(K). Finally, Eq. fl5"2| ) en- 
sures that the Gaussian first term in fl48[) correctly reproduces the variances (x^x^) of 
the original emission function, in particular its r.m.s. widths in the various space-time 
directions. 

Inserting the decomposition (EST) into Eq. (B2|) we obtain for the correlation function 



C(q, K) = 1 + exp[-(f q v {x^x v ){K)\ + SC(q, K) . (53) 

The Gaussian in q results from the Fourier transform of the Gaussian contribution in 
(fj8|); the last term SC receives contributions from the second term SS in P8|) which 
contains information on the third and higher order space-time moments of the emission 
function, like sharp edges, wiggles, secondary peaks, etc. in the source. It is at least 
of fourth order in q, i.e. the second derivative of the full correlator at q = is given 
exactly by the Gaussian in (|53|) . Please note that the exponent of the correlator contains 
no term linear in g; since the correlator must be symmetric under q — > —q because 
it does not matter which of the two particles of the pair receives the label 1 or 2, a 
linear g-dependence could only arise in the form exp(—R\q\). The only type of emission 
function yielding such a g-dependence of the correlator would be a spherically symmetric 
Lorentzian. Any emission function which at large x falls off faster than 1/x 2 results in 
the leading Gaussian behaviour (|53|) instead. This settles, in my opinion, the old issue 
whether Gaussian or exponential fits of the correlation function should be preferred. 

[In the past it has repeatedly been observed that the correlation data appear to be 
better fit by exponentials than by Gaussians. However, as far as I know, this happened 
always when one tried to fit the correlator as a function of the single Lorentz invariant 
variable Qf m = (g ) 2 — q 2 . Contemplating the structure of Eq. (|53|) one realizes that 
such a fit does not make sense: the generic structure of the exponent, —q fJ, q u (x fJj x u ), 
tells us that the term (g ) 2 should come with the time variance of the source while the 
spatial components (g 1 ) 2 should come with the spatial variances of the source. Since all 
variances are positive semidefinite by definition, it does not make sense to parametrize 
the correlation function by a variable in which (g ) 2 and q 2 appear with the opposite 
sign! Such a fit could only work if the time variance and all mixed variances would 
vanish identically, and all three spatial variances were equal. This is certainly not the 
general case in nature. The good exponential fits of the correlation functions from pp 
and e + e~ collisions are thus, in my mind, purely accidental and an empirical curiosity 
without physical meaning. The variable Q irw should not be used for fitting HBT data] 

Please note also that Eq. (|3|) has no factor | in the exponent. If the measured cor- 
relator is fitted by a Gaussian as defined in (|53|), its g-width can be directly interpreted 
in terms of the r.m.s. widths of the source in coordinate space. Any remaining factors 
of \/3, \/3, or \/5 (which you can sometimes find in the literature) are due to reexpress- 
ing the r.m.s. width of the source in terms of certain other width parameters chosen for 



the parametrization of the source in coordinate space. The confusion connected with 
such factors is easily avoided by always expressing the source parametrization directly 
in terms of r.m.s. widths. 



Eqs. (|4q) and fl53|) would, of course, not be useful if the contributions from 5S 
and SC were not somehow small enough to be neglected. This requires a numerical 
investigation. It was shown numerically in Ref. [ |19| that in typical (and even in some 
not so typical) situations SS has a negligible influence on the half width of the correlation 
function. It contributes only weak, essentially unmeasurable structures in C(q,K) at 
large values of q. The reader can easily verify this analytically for an emission function 
with a sharp box profile; the results for the exact correlator and the one resulting 
from the Gaussian approximation (|48 ) are given in [ [18|] and differ by less than 5% in 



the half widths; the exact correlator has, as a function of q, secondary maxima with 
an amplitude below 5% of the value of the correlator at q — 0. We have checked 
that similar statements remain even true for a source with a doughnut structure, i.e. 
with a hole in the middle, which was obtained by rotating the superposition of two 1- 
dimensional Gaussians separated by twice their r.m.s. widths around their center. The 
only situation where these statements require qualification is if the correlator receives 
contributions from the decay of long-lived resonances; this will be discussed in Lecture 
3. 

From Eq. (53) we conclude that the two-particle correlation function measures the 



second order space-time variances of the emission function. That's it - finer features of 
its space-time structure (edges, wiggles, holes) cannot be measured with two-particle 
correlations. The variances (x^x,;) are in general not identical with our naive intu- 
itive notion of the "source radius": unless the source is stationary and has no x-K- 
correlations at all, the variances depend on the momentum K of the pair and cannot 
be interpreted in terms of simple overall source geometry. Their correct physical in- 
terpretation [ |2T], |TJ]] is in terms of "lengths of homogeneity" which give, for each 
pair momentum K, the size of the region around the point of maximal emissivity x(K) 
over which the emission function is sufficiently homogeneous to contribute to the corre- 
lation function. Thus HBT measures "regions of homogeneity" in the source and their 
variation with the momentum of the particle pairs. As we will see, the latter is the key 
to their physical interpretation. 

3.4. Gaussian parametrizations for the correlation function 

A full characterization of the source in terms of its second order space-time vari- 
ances requires knowledge of the 10 parameters (x^Xp). These quantities appear in the 



expression (j42|) for the correlation function but this expression still uses all four com- 
ponents of the relative momnetum g M . However, as already noted only three of the four 
components are independent, due to the mass-shell constraint (|45|). Thus only 6 linear 



combinations of the variances (x ll x u }(K) are actually measurable [[18 



If the source is azimuthally symmetric around the beam axis, this counting changes 
as follows: Even if the source is azimuthally symmetric in coordinate space, the emission 
function S(x, K) in phase space is for finite K no longer azimuthally symmetric because 
the transverse components K± of the pair momentum distinguish a direction transverse 
to the beam direction. There remains, however, a reflection symmetry with respect to 
the plane spanned by K and the beam axis. If we call the direction orthogonal to this 
plane y, all mixed variances which are linear in y must vanish due to this reflection 
symmetry, and the correlator must be symmetric under q y — > —q y . Thus only 7 non- 
vanishing variances (x^x u ) survive in general, of which, due to the mass-shell constraint 



( fi^D only 4 linear combinations are measurable. 

Before the correlator (|53D can be fit to experimental data, the redundant compo- 
nents of q must first be eliminated from the exponent of the Gaussian. At this point 
it is useful to introduce a cartesian coordinate system with z along the beam axis and 
K lying in the x-z-plane. Customarily one labels the ^-component of a 3-vector by 
/ (for longitudinal), the x-component by o (for outward) and the ^/-component by s 
(for sideward). The above choice of the orientation of the x and y axes is natural for 
azimuthally symmetric collision events or event samples because then the transverse 
components K± of the pair momentum distinguish a direction in the transverse plane, 
and it is convenient to orient one of the coordinate axes along this direction such that 
K has only one transverse component: 



K = (K x ,K y ,K z 



(K ± ,0,K L ) 



(54) 



For sources without azimuthal symmetry, e.g. from collisions at finite impact param- 
eter which have been selected according to the orientation of the collision plane, it is 
probably more useful to orient the x axis along the collision plane; then K will be 
characterized by three parameters, Kl, K± and the azimuthal angle $ of K± relative 
to the x-z collision plane. 

A useful formalism for HBT interferometry of finite impact parameter collisions has 
not yet been developed. I will therefore limit my discussion to azimuthally symmetric 
event samples and exploit the symmetry to orient the x-axis along K^. Then from 
(f45|) we see that (3 S = such that 



Q° = f3±q + faqi 



(55) 



with /3j_ = \Kj_\/K° being (approximately) the velocity of the particle pair transverse 
to the beam direction while (3i is its longitudinal component. 

This constraint can now be used in various ways to eliminate the redundant q- 
components from the exponent of Eq. (0). But whichever choice one makes, all the 
.ff -dependent parameters ("HBT radii") in the resulting Gaussian function of q can 
be easily calculated from the variances (x^x"), i.e. by simple quadrature formulae, for 
arbitrary emission functions S(x,K). The relation between the HBT parameters and 
the variances is model-independent, i.e. it does not depend on the form of the emission 
function S(x, K). 

Here I discuss two specific parametrizations: the standard Cartesian one (mostly 
iiii 



for historic reasons 
Podgoretskh one [ p|, H, |24 



and the more physically motivated Yano-Koonin- 



3.4-1- Standard Cartesian parametrization 



The standard form 
obtained by using 



14], ^(J for the parametrization of the correlation function is 
) to eliminate q° from Eq. 



C(q, K) = 1 + exp 



One obtains 



£ R%{K)q iqj 

i,j=s,o,l 



(56) 



where the 6 HBT radius parameters Rij are defined in terms of the following variances 



of the source function [ 14, EQ, 28 1: 



R l( K ) = ((£< - Pit)(xj ~ Pjt)) , i,j = s, o, I 



(57) 



For an azimuthally symmetric sample of collision events C(q,K) is symmetric with 
respect to q s — > — q s [ pT8[] . Then R 2 S = R 2 sl = and 

C(q, K) = l + exp [-R 2 s (K)q 2 s - R 2 (K)q 2 - R 2 (K)q 2 - 2R 2 ol (K)q qi\ , (58) 

with 

(59) 
(60) 
(61) 
(62) 

The cross-term (|B~2"D was only recently discovered [ p0[| . Clearly these HBT radius 
parameters mix spatial and temporal information on the source in a non-trivial way. 
Furthermore, since they multiply combinations of the components q^ which are not 
invariant under longitudinal boosts of the measurement frame, their interpretation 
depends on the frame in which the particle momenta are specified. This complicates 
their physical interpretation. An extensive discussion of these parameters, in particular 
of the meaning of the generally non- vanishing cross-term R 2 ol , can be found in Refs. [ 
O, [TR [TP], 120, 123, RTJ, where the expressions (153)- (K21) were analyzed analytically and 
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Rf 
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= ((2- 


Pit) 2 ) , 


Rl 


:k) 


= ((*- 


P±i)(z - Pit)) 



numerically for a large class of (azimuthally symmetric) model source functions. Some 
of these results will be discussed in the Lecture 3. 

An important observation resulting from these studies is that the difference 

R 2 iS = R 2 Q -R 2 = P 2 ± (t 2 ) - 2P ± (xt) + ((x 2 ) - (y 2 )) (63) 

is generally dominated by the first term on the r.h.s. and thus provides access to the 



lifetime At = J (t 2 ) — (t) 2 of the source (more exactly: the duration of the particle 
emission process) [ |3T[ . However, in relativistic heavy- ion collisions, due to rapid ex- 
pansion of the source one would not expect (t 2 ) to be generically much larger than 
either (x 2 ) or (y 2 ) unless there is a phase transition to a quark-gluon plasma and the 
collision fireball is initiated within a certain range of energy densities above the critical 



energy density where the transition occurs [32]. In the standard fit one is not sensitive 



to small values of At since Eq. (p3|) then involves a small difference of two large num- 
bers, each associated with standard experimental errors. The factor /3j_ < 1 in front of 
{t 2 ) further complicates its extraction, in particular at low K± where At(K) is usually 
largest (see below). Successful attempts to determine the duration of particle emission 
from HBT measurements have been reported from low-energy heavy-ion collisions (us- 
ing 2-proton correlations) where the measured lifetimes are very long: 25 ± 15 fm/c in 
Ar+Sc collisions at E/A = 80 MeV [ and 1400 ± 300 fm/c in Xe+Al collisions at 
E/A — 31 MeV [ |34J (the latter is the typical evaporation time of a compound nucleus). 
Two-pion correlations at ultra-relativistic energies (E/A = 200 GeV) so far failed to 
yield posit iveevidence for a non- vanishing emission duration [ |35|, |36|] , except for the 



heaviest collision system Pb+Pb [ 3"7| , but even there the effective lifetime is only a few 
fm/c. 

3.4-2. Yano-Koonin-Podgoretskii (YKP) parametrization 



The Yano-Koonin-Podgoretskh parametrization of the correlation function is the 
generalization to azimuthally symmetric systems of the Yano-Koonin parametrization 
[ It was first written down by M.I. Podgorestskh in 1983 for moving azimuthally 



symmetric, but non-expanding sources [ [39]], with iC-independent parameters. Not 



knowing about this paper, we reinvented it in [ 18f , with independent parameters for 



expanding azimuthally symmetric systems. The YKP parametrization is based on an 



elimination in Eq. (|53| ) of q and q s in terms of q± = y q 2 + q 2 , q°, and g 3 , using Eq. (|55|): 
C(q, K) = 1 + exp [-it!! Ql ~ R\ U ~ i<f?) ~ fao + ^J) (? • • (64) 



Like the standard Cartesian form (^) it has four iT-dependent fit parameters, but now 
only three of them, R±(K), Rn(K), and Ro(K), have dimensions of length while the 
fourth parameter, U(K), is a 4- velocity with only a longitudinal spatial component: 

U(K)=j(K)(l,0,0,v(K)), with 7 = ^=i==. (65) 



This parametrization has the advantage that the "YKP radii" R±, R\\, and Rq extracted 
from such a fit do not depend on the longitudinal velocity of the observer system in 
which the correlation function is measured; they are invariant under longitudinal boosts. 
Their physical interpretation is easiest in terms of coordinates measured in the frame 
where v(K) vanishes. There they are given by [ |i~8| 



R\{K) = R 2 s (K) = (y 2 }, (66) 
R\{K) = ((z-Ax//3 ± ) 2 )-A 2 (y 2 )//5i^(5 2 ), (67) 

R%(K) = ^(t-x/(3 ± ) 2 y(f}/(3l^(P} } (68) 
where in the last two expressions the approximation consists of dropping terms which 



were found in [ [I8j to be generically small. (A more quantitative discussion of this point 
will follow in Lecture 3.) The first expression fl6"6] ) remains true in any longitudinally 
boosted frame. 

Eq. (pH) shows that the YKP parameter Rq(K) essentially measures the time 



duration At(K) = J (t 2 ) during which particles of momentum K are emitted, in the 
frame were the YKP velocity v(K) = 0. It enters as the leading contribution in Rq, is 
fitted directly and no longer obtained as the difference of two large fit parameters as in 
the standard Cartesian fit. 

Eqs. (B^)-(B^) were written down in the special frame where v(K) = which 
we call the Yano-Koonin (YK) frame. In Refs. [ 0, |24[] it is shown that for a large 
class of models this frame essentially coincides with the longitudinal rest frame of the 
fluid cell around the point x(K) of maximum emissivity at momentum K (i.e. the 
Longitudinal Saddle Point System LSPS [ P2|Q. This is true also for sources which 
are not longitudinally boost-invariant and for which the LSPS and the LCMS (the 
Longitudinally CoMoving System in which the pion pair has f3i = [ [H|]) do not 
coincide. 

In general the measurement system will not coincide with the (iT-dependent) YK- 
frame, and the YKP radii will be given by more complicated combinations of the 
space-time variances of the source expressed in the coordinates of the measurement 
frame. This is the simple result of a Lorentz-boost between the two frames. However, 
I stress that in any frame the YKP parameters can again be easily calculated from 
the second order moments of the source function S(x,K), i.e. by simple quadrature. 
Introducing the notational shorthands 



(69) 



B = ((i-Ae/^) 2 ) , (70) 

C = ((t-e//3±)(«-Af//?±)), (71) 

where £ = x + and (y) = (xy) = for azimuthally symmetric sources such that 
(£ 2 ) = (x 2 ~ y 2 ), one finds in an arbitrary longitudinal reference frame 



A + B 




(72) 



R 2 = - ( ^(A + E) 2 -^ 2 - A + = B-uC, (73) 
1 



R 2 = ^U{A + Bf-AC 2 + A- Bj = A-vC . (74) 

The Yano-Koonin velocity v is zero in the frame where the expression flTT| ) for C vanishes 
[ [L8|]; this fixes also the sign in front of the square root in (0). For small values of C 
the Yano-Koonin velocity is given approximately by 

C 

(75) 



A + B (t 2 ) + (P) ' 



where in the second approximation we again neglected generically small terms [ 18 



proportional to (zx), (xi), and {x 2 — y 2 ). The accuracy of the approximate expression 
( [75]) for v(K) was tested numerically in [ 24 1 and found to be excellent in the situations 
discussed below. In the same limit the expressions for Rq and R 2 simplify to R 2 , ~ A 
and R 2 w B, in agreement with (67) and (|6q) . 

Since the standard Cartesian and YKP parametrizations are mathematically equiv- 
alent (being simply based on a different choice of independent g-components) , the re- 
sulting HBT parameters obey simple relations [ p3|j : 

R 2 S = R±, (76) 

Rle = Rl' R 2 = {Rl + v 2 R\) , (77) 

R 2 = (l-(3f)Rl + j 2 (p l -vf(Rl + R^ J (78) 

R 2 0l = P ± (-Ai?| + 7 2 (A - (^o + ^)) • (79) 

Although a mathematical triviality, these relations provide a powerful consistency 
check on the experimental fitting procedure, of similar value as the relation [ [TB|, |H| 
UmK ± ^o{R {K) — R S (K)) = which results from azimuthal symmetry. 



4. LECTURE 3: HBT CORRELATIONS FOR EXPANDING SOURCES 

In this third lecture I will present a quantitative discussion of the HBT parameters, 
both in the standard Cartesian and in the YKP fits. The emphasis will be on their 
^-dependence because, as discussed in the previous two lectures, only a careful study 
of this pair-momentum dependence permits a separation of the geometric from the 
dynamic aspects of the source. You probably remember that at the beginning of Lecture 
2 I stressed that a model-independent reconstruction of the emission function from HBT 
measurements is not possible. Therefore any quantitative discussion must necessarily 
occur within the framework of specific source models. 



I will choose a relatively simple analytical parametrization of the emission function 
which was first suggested by T. Csdrgo and B. Lorstad in 1994 in an unpublished 
preprint (see also [ ^2], [2(], [L4|, |IB|| ) and which incorporates many of the (as we believe) 
relevant physical features of the typical sources created in relativistic nuclear collisions: 
approximate thermalization at decoupling, finite transverse and longitudinal extension 
of the source, collective expansion in both the longitudinal and transverse directions, 
and a finite duration of the particle emission process. All these features are controlled 
by parameters which can be freely tuned, thus allowing for extensive parameter studies 



19 



, |30| which have given us a good intuitive understanding of which properties 
of the source are important for the correlation function and which aren't, and how to 
look in the correlation function for specific space-time properties of the emitter. 



4.1. A simple model for an expanding source 



Let us consider the following model for the emission function [ IS, 22 



with 



and 



S(x,K) 



H(x) 



2 J 



(2vr) : 



K ■ n(x) exp 



K ■ u(x) — fi(x) 



71 At 



cxp 



T(x) 

(V - Vo) 2 (r - Tof 



H(x) 



2R 2 2(Ar]y 



2(At)< 



(82) 



K ■ n(x) = Mj_ cosh(?7 — Y) . 

The factor 2 J + 1 counts the spin degeneracy of the observed particle species and is 
included because the detectors usually do not identify the polarization of the measured 
particles. There is no such factor for isospin because of the requirement that the two 
particles in the pair be identical, i.e. have the, for example, the same electric charge. 
The Lorentz covariant Boltzmann factor exp[— (K-u(x) — fi(x))/T(x)] implements the 
assumption of local thermal equilibrium of the emitted particles at freeze-out, with 
local temperature T(x) and chemical potential n(x), and the collective expansion of 



the source with hydrodynamic flow 4-velocity 



[x 



I will here take T and /i as 



constants and refer to Refs. [ |2^, Q for an investigation of the effects of temperature 
gradients. The factor H(x) describes the geometric properties of the source; it can be 
interpreted as a space-time modulation of the fugacity exp[/i(x)/T(x)]. Space-time is 
parametrized by longitudinal proper time r = \Jt 2 — z 2 and space-time rapidity r\ = 
\ ln[(t + z)/(t — z)] in the longitudinal and temporal directions, and by r = ^/x 1 + y 2 
and <p in the transverse directions. The volume element is then d 4 x = r dr drjr dr d(f). 
The Gaussian factors exp[— r 2 /(2R 2 )] and exp[— (r/ — i] ) 2 / (2(Arj) 2 )] provide smooth 
geometric limits for the source in the transverse and longitudinal directions, scaled by 
R and L = rArj, respectively. The function H(x) is normalized to the total comoving 
3- volume according to 



d A xH(x) 



nr r 



where 



2R 2 



+ y; 



2r r/ r 



At] 



(83) 



34) 



are the r.m.s. expectation values for the transverse radius r and for 77, respectively. 
(The r.m.s. widths for x and y are each given by R.) If the Gaussians in H(x) were 
replaced by box functions, the equivalent box dimensions (with the same r.m.s. radii) 
would be R = 2R, fj = \^3Ar]. (Here you see the "unneccesary" factors \[2 and \/3 
mentioned before!) 



K-n(x) is the flux factor through the freeze-out hypersurfaces whose normal di- 
rection is given by the unit vector In our model these hypersurfaces are for 
simplicity assumed to be surfaces of constant longitudinal proper time r, parametrized 
by surface coordinates E( T )(x) = (r cosh 77, r cos 0, rsin0, Tsmhrj). The last factor in 
H{x) provides a Gaussian smearing of the proper time around the mean values tq with 
dispersion Ar, thereby implementing particle emission over a finite time interval of 
order Ar in the local comoving frame. With these assumptions the flux factor reduces 
[ 18 1 to the form given in Eq. (|S2D. For Ar — > the Gaussian in r approaches a 8- 



function centered at To, simulating instantaneous freeze-out at constant proper time as 
often implemented in hydrodynamical situations. 

[Let me make some side remarks on hydrodynamical simulations here, because they 
are a very popular tool for computing the space-time dynamics of the reaction zone 
in heavy-ion collisions. In hydrodynamics freeze-out is always assumed to occur on a 
sharp hypersurface (not a smeared one as here), and one writes the emission function 
in the form [ J5| 

= 2J+1 f K^dh,{x')5^{x-x') 
K ' J (2tt) 3 k exp[{K-u{x') - fi{x'))/T{x')} ± 1 ' 1 ; 

where d^a^ix') is the normal vector on the freeze-out surface S(x'), and we have cor- 
rectly accounted for the quantum statistical effects in the local thermal distribution 
functions. The latter are unimportant for heavy particles but should, in a quantitative 
comparison with data, be included for pions - at least in the single-particle spectrum. 
Inserting the ansatz (|85|) into the expression (^3) for the single-particle spectrum one 
obtains the well-known Cooper-Frye formula [ |4lj 



d 3 K Jx 

with the local equilibrium distribution function 



dN 

E k^77= I K»d A a,{x)f{x,p) (86) 



(2vr) 3 exp[(K-u(x) - fi(x))/T{x)} ± 1 ' ^ 



For the numerator of the correlator in fl39|) one obtains 

K»d 3 a^x) K v d\ v {y) f(x, K) f(y, K) exp[zg-(a; - y)\ 



s 



similar, but not identical with to the one given in [ Q . In [ Q each of the two distri- 
bution functions under the integral featured on-shell arguments p a and p&, respectively, 
instead of the common (off-shell) average argument K as in (|55|). This error in [^| can 
be traced back to an inaccurate transition from finite discrete volumes along the freeze- 
out surface S to the continuum limit [|43|]. Taking over this inaccuracy produces (in 
particular for very rapidly expanding sources) unphysical oscillations of the correlation 
function around unity at large values of q (see e.g. [ which are inconsistent with the 
manifestly positive definite nature of the exchange term noted in Lecture 1. Recently 
Aichelin [ [fij pointed out that, for technical reasons, the same erroneous assumption is 
made in all simulations of HBT correlations using Monte Carlo event enerators. This 
problem should certainly receive more attention in the future.] 

Since according to Eq. ( pETp the time-component of the pair momentum may be 
approximated by the on-shell value Kq = Ek = \/m? + K 2 , the pair momentum K can 



be parametrised in terms of the momentum rapidity Y = \ ln[(l + /3/)/(l — A)] an d the 



transverse mass M± = ym 2 + Kj_, 

K» = (M ± cosh Y, K ± , 0, M ± sinh F). (89) 

We implement longitudinal and azimuthally symmetric transverse expansion of the 
source by parametrising the flow velocity in the form 

/ x y \ 

u^(x) = cosh rji (r, if) cosh rj t (r), — sinh?7t(r), — sinh r]Ar), sinh rji (r, rj) cosh r]Ar) ) . 
V ' t r J 

(90) 

For the longitudinal flow rapidity we take rji (r, rj) — r\ independent of r, i.e. we assume 
a Bjorken scaling profile [ 45] Vi — z/t in the longitudinal direction. The growth of the 



longitudinal flow velocity in the longitudinal direction is then automatically limited by 
the Gaussian in rj in (|31~D. For the transverse flow rapidity we take a linear profile of 
strength r\f. 

Vt(r)=Vf(^)- (91) 

The scalar product in the exponent of the Boltzmann factor can then be written as 

x 

K ■ u(x) = M_|_ cosh(r/ — Y) cosh T] t (r) — Kj_— sinh^ t (r) , (92) 

Please note that for non-zero transverse momentum K±, a finite transverse flow breaks 
the azimuthal symmetry of the emission function via the second term in fl92"|). For rjf = 
the source has no explicit K±- dependence, and M±_ is the only relevant scale. As will 
be discussed later this gives rise to perfect M^-scaling of the YKP radius parameters 
in the absence of transverse flow, which is again broken for non-zero transverse flow [ 

g. 

Besides rjf, the model parameters are the freeze-out temperature T, the transverse 
geometric (Gaussian) radius R, the average freeze-out proper time r as well as the 
mean proper emission duration At, the centre of the source rapidity distribution rjo, 
and the (Gaussian) width of the space-time rapidity profile Ar]. A rough spatial picture 
of the source at various fixed coordinate times can be gleaned from Figs. 1 and 2 in 



Ref. [ f|7| to which I would like to refer those readers having trouble with visualizing 



Gaussians in rj and r in regular Cartesian coordinates. Although the source in Ref. 



47j has sharp edges in space and time whereas ours is smoothed by Gaussian profiles, 
the essential space-time features of the sources are very similar. 

The calculations presented below were done for pions (m = m n ± = 139 MeV/c 2 ) 
and kaons (m = m^± = 494 MeV/c 2 ) and (except were noted otherwise) for the fixed 
set of parameters T = 140 MeV, R = 3 fm, r = 3 fm/c, Ar = 1 fm/c, A77 = 1.2, 
and r] = (i.e. all our rapidities will be given relative to the rapidity i] of the cm. 
of the source). The strength rjf of the transverse flow will be varied systematically 
to investigate its effects on the HBT parameters. A detailed discussion of how the 
variation of some of the other parameters affects the correlation function can be found 



in Refs. [ 24, 30 



4.2. independence of the HBT radii 

4-2.1. Analytical approximations - HBT radii as lengths of homogeneity 



For a qualitative understanding of the physical origin of the pair-momentum de- 
pendence of the HBT parameters it is instructive to start from their model-independent 



expressions in terms of space-time variances and evaluate the latter approximately ana- 
lytically. For the standard Cartesian fit parameters ( |59|) -(l6^) this was done, at different 



levels of accuracy, in Refs. [ [T3|, [19|, |2T], by exploiting the method of saddle-point 
integration for the evaluation of the variances. This method was introduced by Makhlin 



and Sinyukov [ 25] in the context of infinitely long sources with boost-invariant longi- 
tudinal expansion where they used it to show that the longitudinal HBT radius Ri is 
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Fig.l. i^-dependence of the standard Cartesian HBT radii R Q (top), R s (middle), 
and Ri, evaluated in the LCMS (the frame where Kj, — 0), for an infinitely long 
source (A77 — > 00 in Eqs. ( |8C)| , |Sl| ) ) with boost-invariant longitudinal expansion and a 
linear transverse flow rapidity profile. The strength of the tranverse flow rjf is varied 
between rjf = and rjf = 0.9 as indicated. The duration parameter At was set to 
(locally instantaneous freeze-out). For such a source the cross-term R a i vanishes 
identically in the LCMS. Solid lines represent results from a full numerical evaluation 
of the integrals in Eqs. (|59|-|62|). Dashed lines give the results from the leading-order 
saddle-point approximation to these integrals. One sees that the saddle-point integra- 
tion misses the rise of R at small K± (i.e. the contribution from the finite duration 
of particle emission in the LCMS frame, see text) and gives a very bad approxima- 
tion to Ri at small K±. The solid line for Ri for 77/ = is given analytically by 
Rl = t y/T/M± sjK 2 (M± /T) /Ki (M± /T) [|| while the saddle-point approxima- 
tion (dashed) gives for rjf — the Makhlin-Sinyukov formula Ri = tq^JTJM±_ [ |25| . 
(Figure taken from Ref. [ IHj.) 

finite and determined by the inverse of the longitudinal velocity gradient, noting for 
the first time that Ri has the property of a "longitudinal length of homogeneity" in the 
source rather than being related to the longitudinal geometric size of the entire source. 



It was later observed [ |14], |i| that not all of the important qualitative features of the 
independence of the HBT parameters can be obtained from the leading term in the 
saddle point approximation (see Fig. 1). In particular, in the presence of transverse 
flow the saddle point moves away from the beam axis r = 0, and this must be taken into 



account in order to obtain reasonable approximations [ [19] . Unfortunately, this renders 



the whole procedure rather cumbersome [ [L9|, and in the end, e.g. for a quantitative 



comparison with data, a full numerical evaluation of the integrals for the variances 
cannot be avoided. 

In spite of their unreliability on a quantitative level, the analytical results from 
saddle-point integration are still very instructive on a qualitative level. I will here 
discuss the leading results for a longitudinally finite (L = r Ar]) source of type QgDD in 
the limit At = 0, in the LCMS (i.e. in the Longitudinally CoMoving System where 
the pairs have Kl = 0). One finds [ [H] 

R 2 s (M ± ,K L = 0) = Rl, (93) 
R 2 o (M ± ,K L = 0) = Rl + Pl(AQ 2 , (94) 
R 2 (M ± ,K L = 0) = Ll, (95) 

R 2 ol (M ± , K L = 0) = AU , (96) 



(97) 
(98) 
(99) 



are functions of M± 


defined by 


1 


1 1 

R 2 + Rj 


m~ 


1 


1 1 




AU = V2 Ur 2 - 


f L 2 ~r ) 



with the transverse and longitudinal "dynamical lengths of homogeneity" 



R «^=^Jwridwr'lk- (100) 



M mj = To ^L = _L^L (ion 



Strictly speaking, expression ( |100| ) is only valid for weak transverse flow r/f <^ I. In 
( |101|) we have defined the longitudinal flow 4-velocity u\ = (cosh^, 0, 0, sinhr^) = 
(cosh 77, 0, 0, sinh?]). 

The physical interpretation of these results is quite interesting: in addition to ge- 
ometry (implemented by the Gaussian cutoff factors in the function H(x)), dynamics 
affects the HBT radii through the dynamical homogeneity lenghts Rh, Lh- The latter 
are inversely proportional to the gradients of the expansion velocity field in the respec- 



tive direction, but smeared by a thermal smearing factor yT/M± resulting from the 
random thermal motion of the particles around the fluid velocity. The HBT radii are 
determined by the shorter of the two (geometric or dynamic) lengths scales. In the 
absence of random thermal motion (T — > 0) any velocity gradient in the system would 
lead to a vanishing dynamical length of homogeneity and consequently to vanishing 
HBT radii. At finite T, the dynamical smearing decreases with increasing transverse 
mass Mj_, leading to a decrease of the HBT radii at large M±. (It turns out that the 
-scaling of the HBT radii at large Mj_ suggested by these analytical expressions 



sJl/M x 



is unreliable and a consequence of the saddle-point approximation [ [L9[. A numerical 



analysis [ |46|, |J] shows that the power of M± itself, by which the HBT radii decrease 
for increasing M±, is proportional to the expansion velocity gradient.) 

The following intuitive picture results from these considerations: if the expansion 
velocity is small, i.e. all velocity gradients can be essentially neglected over the range 
where the geometric Gaussians in H(x) are large, then HBT measures the geometric 
parameters R, L which tell you where the function H(x) (and thus the whole emission 
function S(x,K)) gets cut off. If, on the other hand, the velocity gradients are large, 
they effectively cut off the emission function at a distance Rh resp. L H from the saddle 
point, and the matter outside these homogeneity regions decouples from the correlator 
because it cannot contribute particles with sufficiently small relative momenta to see 
the effects of quantum statistics. This explains my statement above that HBT does in 
general not measure the geometry of the source, but rather the regions of homogeneity 
inside the source at a given wavelength 1/K. 

Of course, a space-time dependence of the temperature field T(x) can induce ad- 
ditional gradients into the emission function and thus affect the size of the regions of 
homogeneity in the source and the HBT radii. This was investigated in some detail in 
Refs. f 1211 RDI to which I refer the interested reader. 



The last point to be discussed is the origin of the quantity At* in Eq. ([J|) for 
R 2 . Comparing with Eq. fl63|) we see that it has the meaning of an effective source 
lifetime. But where does it come from, since we set the width At of the proper time 
distribution to zero? This apparent paradox has an interesting answer which is reflected 
in the mathematical structure of Eq. (^): since the correlator receives non- vanishing 
contributions from a longitudinal region of homogeneity of size Ri = L*, it probes 



emission from different points z at different times t = yr^+z 2 along the proper-time 
hyperbola r = To, with maximal range — <z<L*. Thus even for sharp freeze-out at 
constant proper time the correlator sees a non-vanishing effective lifetime in the fixed 
observer frame (here the LCMS) which is in principle measurable via the difference 
R 2 Q — R 2 S . Since L* is a decreasing function of M±, so is At*, and for large Mj_ the 
difference R 2 — R 2 vanishes. (If At were nonzero, the difference would at large M± 
approach the limit (At) 2 .) 



4-2.2. Cartesian HBT radii in the CMS and LCMS 



In Fig. 2 I show the HBT radius parameters from the standard Cartesian fit ([58]) 
for pion pairs with cm. rapidity Y = 1.5 where the fit of the correlator is done in 
the CMS [ 0. The different thick curves correspond to different strengths rjf of the 
transverse flow. Without transverse flow R s is Mj_-independent because the we consider 
(jS0|) for constant temperature and neglect possible transverse temperature gradients. 
As the transverse flow increases, R s develops an increasing dependence on M±. As 
shown in Fig. 6 below it can be approximated by an inverse power law, with the power 
increasing monotonously with rjf [ [19], g4]| . 

Ri features a very strong Mj_-dependence even without transverse flow, due to the 
strong longitudinal expansion of the source. It can also be described by an inverse 
power law, with a larger power ~ 0.55, in rough agreement with the approximate 
^/r/Mi-scaling law suggested in [ [25[ (see, however, [ [19], |28j for a more quantitative 
discussion). The increase of R a at small Mj_ is due to the contribution fl63|) from the 
effective lifetime. As seen in Fig. 5 below, in the YK frame (source rest frame) the 
latter is of order 2.5 fm/c at small M^; Fig. 2b shows that its effect on R compared to 
R s in the CMS is much smaller (and thus more difficult to measure). Fig. 2d shows that 
the cross-term is small in the CMS but non-zero. It vanishes at K± = by symmetry 
and also becomes small again at large Kj_. 
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Fig. 2. The standard Cartesian parameters R s (a), R Q (b), Ri (c) 



WUMeV) 



and R ol (d) in the 



CMS for pion pairs with cm. rapidity Y = 1.5, as functions of M± for 3 different 
values for the transverse flow rjf. The thick lines are exact numerical results from 
Eqs. (|59|-|62|) , the thin lines are obtained from the analytical approximations given in 
Ref. (Figure taken from Ref. [|3(J.) 

The thin lines in Fig. 2 show for comparison approximate results for the HBT radii 
calculated from the approximate analytical results given in Ref. [ Q which were derived 

j2D by saddle point integration. It is clear that this method fails 
for a quantitative discussion of this approximation), and that the 



by evaluating Eqs. (|59 
here 



see Ref. [ 

analytical expressions should not be used for a quantitative analysis of HBT data. 

Fig. 3 shows the same situation as Fig. 2, but now all HBT radii are evaluated in 
the LCMS (longitudinally comoving system [ 3l| ) which moves with the pair rapidity 



Y = 1.5 relative to the CMS. A comparison with Fig. 2 shows the strong reference 
frame dependence of the standard HBT radii. In particular, the cross-term changes 
sign and is now much larger. The analytical approximations from Ref. [ |22| work much 
better in the LCMS [ g2|, but for R Q and R 2 ol they are still not accurate enough (in 
particular in view of the delicate nature of the lifetime effects on R a ). 

4-2.3. The Yano-Koonin velocity 



Fig. 4 shows (for pion pairs) the dependence of the YK velocity on the pair mo- 
mentum K. In Fig. 4a we show the YK rapidity Y YK = ^ln[(l + v)/(l — v)] as a 
function of the pair rapidity Y (both relative to the CMS) for different values of K±, in 
Fig. 4b the same quantity as a function of K± for different Y. Solid lines are without 
transverse flow, dashed lines are for rjf = 0.6. For large K± pairs, the YK rest frame 
approaches the LCMS (which moves with the pair rapidity Y); in this limit all pairs 
are thus emitted from a small region in the source which moves with the same longi- 
tudinal velocity as the pair. For small K± the YK frame is considerably slower than 
the LCMS; this is due to the thermal smearing of the particle velocities in our source 
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Fig. 3. Same as Fig. 2, but now evaluated in the LCMS. Please note the change of 
sign and magnitude of the cross- term. (Figure taken from Ref. [ p0|.) 

around the local fluid velocity u^ix) [ The linear relationship between the rapidity 
Y YK of the Yano-Koonin frame and the pion pair rapidity Y is a direct reflection of 
the boost- invariant longitudinal expansion flow [ For a non-expanding source Y YK 
would be independent of Y. Additional transverse flow is seen to have nearly no 
effect. The dependence of the YK velocity on the pair rapidity thus measures directly 
the longitudinal expansion of the source and cleanly separates it from its transverse 
dynamics. A detailed discussion of these features is given in Ref. [ [?3J where it is also 
shown that the YK velocity is always very close to the velocity of the Longitudinal 
Saddle Point System LSPS (i.e. to the longitudinal velocity of the fluid element around 
the point of maximal emissivity at momentum K). This last observation establishes 
the usefulness of the YK velocity (which can be directly extracted from an YKP fit to 



the data) as a measure for the longitudinal expansion velocity of the source. 




Fig. 4. (a) The Yano-Koonin rapidity for pion pairs, as a function of the pair cm. 
rapidity Y, for various values of K±_ and two values for the transverse flow r\j. (b) 
The same, but plotted against K± for various values of Y and rjf. (Figure taken from 
Ref. [HO 

4-2.4- YKP radii: M^-scaling and transverse flow 

In the absence of transverse flow, a thermal source like fl80|) depends on the particle 
rest mass and on the transverse momentum K± only through the combination M]_ = 
m 2 + K\ (see Eq. (|92|)). Furthermore, the source is then azimuthally and x — > —x 
reflection symmetric. Hence (xt), (xz), and (x 2 — y 2 ) all vanish and the approximations 
in Eqs. ([H],^) become exact. As a result, all three YKP radii (|d^)-(|H5|) are only 
functions of M±, too (as well as of Y, of course), i.e. they do not depend explicitly on 



the particle rest mass. 
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Fig. 5. The YKP radii R±, Ru, and Ro (from top to bottom) for vanishing transverse 
flow (left column) and for rjf = 0.6 (right column), as functions of M± for pairs at 
Y cm = 0. Solid (dashed) lines are for pions (kaons). The breaking of the M^-scaling 
by transverse flow is obvious in the right column. Also, as shown in the lower right 
panel, for nonzero transverse flow Rq does not agree exactly with the effective source 
lifetime \J (t 2 ). (Figure taken from Ref. [ p4|.) 

This is seen in the left column of Fig. 5 where the three YKP radii are plotted for 
y cm = pion and kaon pairs as functions of M±; they agree perfectly. The transverse 
radius here shows no M^-dependence due to the absence of transverse temperature 
gradients, but even with temperature gradients it would only depend on M±. (Of 
course, this discussion neglects resonance decays which will be studied in Sec. 4.3.) The 
very strong M^-dependence of the longitudinal radius parameter R\\ is again due to the 
strong longitudinal expansion of the source. Note that Mj_-scaling in the absence of 
transverse flow applies only to the YKP radius parameters: since the expressions (j60|)- 
(|62|) involve nonvanishing variances with f3±- or /3rprefactors (which depend explicitly 
on the rest mass), the HBT radii from the standard Cartesian fit do not exhibit Mi- 
scaling. 

For non-zero transverse flow rjf 7^ this Mj_-scaling is broken by two effects: first, 
the second term in (|92[) destroys the M^-scaling of the emission function itself, and 
second the /3-dependent correction terms in (|67|,58) are now non-zero because the same 
term also breaks, for K± 7^ 0, the x — > —x and x — * y symmetries. The magnitude 
of the associated scale breaking due to the pion-kaon mass difference is seen in the 
right column of Fig. 5 for 77/ = 0.6. The effects are small and require very accurate 
experiments for their detection. However, the sign of the effect is opposite for R\\ and 



for R_i, R Q which may help to distinguish flow- induced effects from resonance decay 
contributions. 

Since for Y cm = the YK and CMS frames coincide, fii = in the YK frame and 
the approximation in fl6~7|) remains exact even for non-zero transverse flow. The same 
is not true for the approximation in (|68|), and therefore I show in the lower right panel 



of Fig. 5 also the effective source lifetime J (t 2 ) for comparison. The apparently rather 
large discrepancies between the YKP parameter R and the effective source lifetime is 
due to a rather extreme choice of parameters: a large flow transverse flow and a small 



intrinsic source lifetime of At = 1 fm/c in fl8T|). Since y (t 2 ) approaches At in the 
limit of large M_l while the dominant [ |24| correction term (x 2 — y 2 ) does not depend 
on At, the YKP parameter R will track the effective source lifetime more accurately 
for larger values of At (and for smaller values of rjf). 



Why do y (t 2 ) and Rq increase at small M±? Due to the rapid longitudinal ex- 
pansion, the longitudinal region of homogeneity R» is a decreasing function Mj_. Since 
for different pair momenta Rq measures the source lifetime in different YK reference 
frames, the freeze-out "hypersurface" will in general appear to have different shapes 
for pairs with different momenta. Only in our model, where freeze-out occurs at fixed 
proper time r (up to a Gaussian smearing with width At), is it frame-independent. It 
is thus generally unavoidable (and here, of course, true in any frame) that freeze-out at 
different points z in the source will occur at different times t in the YK frame. Since a 
^-region of size R\\ contributes to the correlation function, Ru determines how large a 
domain of this freeze-out surface (and thus how large an interval of freeze-out times in 
the YK frame) is sampled by the correlator. This interval of freeze-out times combines 
with the intrinsic Gaussian width At to yield the total effective duration of particle 
emission. It will be largest at small pair momenta where the homogeneity region R» is 
biggest, and will reduce to just the variance of the Gaussian proper time distribution 
at large pair momenta where the longitudinal (and transverse) homogeneity regions 
shrink to zero. The rise of At(K) at small K is thus generic. 

While the strong Mj_-dependence of the longitudinal radius parameter Ru arises 
from the strong longitudinal expansion, the weaker Mi-dependence of the transverse 
radius parameter reflects the weaker transverse expansion of our source. Following a 



suggestion by Th. Alber [ [48|, this relation can be made quantitative: in Fig. 6 we 
plot in the left column the transverse and longitudinal YKP radii R± and Ru versus 
Mi on a double-logarithmic scale. We see that both can be approximately represented 
by power laws. (The same is not true for Rq.) While in such a plot the slope of R± 
clearly increases with the strength r)f of the transverse flow, the slope of R\\ appears to 
be insensitive to transverse flow. This can be seen quantitatively in the right column 
of Fig. 6 where we plot the powers a±,a\\ extracted from a fit 

R±(M ± ) oc M~ a± - , R\\(M ± ) oc M~ a " (102) 

as a function of rjf. As indicated in Fig. 6b the extracted power a±_ for R± depends 
somewhat on the fit region because R± doesn't follow an exact power law; indepen- 
dent of the fit region it increases, however, monotonously and nearly linearly with the 
strength rjf of the transverse flow. Kaons "feel" the transverse flow more strongly than 
pions, as reflected by the somewhat larger powers a± at fixed rjf. Note that even for 
a rather strong transverse flow rjf = 0.6 (heavy-ion data seem to require less flow) a± 



remains below 0.25. 




Fig. 6. (a) R± as a function of M± at Y CM = 0, for pions (solid) and kaons (dashed) 
and different transverse flow rapidities rjf. (b) The scaling coefficient a± defined by 
R± w M ^ a± for pions (solid) and kaons (dashed) as a function of the transverse flow 
rapidity rjt. The different results obtained by fitting in the regions M± — m < 500 
MeV/c 2 and M± — m < 1000 MeV/c 2 are shown separately, (c) Same as (a), but for 
Rh. (d) Same as (b), but for aii. (Figure taken from Ref. [ B.) 

The power ot\\ for R\\ is shown in Fig. 6d. It has a much larger value of 0.55-0.56 
even for rjf — 0, reflecting the strong boost-invariant longitudinal expansion. As rjf 
increases, a\\ also increases, but only by a few percent. (Note the suppressed zero in 
Fig. 6d!) Again, kaons are affected more strongly by the transverse flow than pions, 
but altogether the Mj_-dependence of R» is rather insensitive to transverse flow. 



As observed by Th. Alber in his thesis [ |48| these features agree qualitatively with 



the heavy-ion data: for ot\\ he found values of order 0.4-0.5, while a± remained smaller, 
of order 0.1-0.2; in 32 S-induced collisions, the central values for both numbers showed a 
systematic tendency to increase with the mass of the target nucleus, indicating stronger 
collective flow in larger collision systems than in smaller ones. The error bars on the en's 
were, however, large, and one should wait for independent confirmation before firmly 
drawing such a conclusion. 

The F-dependence of the Yano-Koonin rapidity Y YK and the M^-dependence of 
Rn can thus be used as a quantitative measure for the longitudinal expansion of the 
source which is hardly affected at all by the presence and strength of transverse flow. 
On the other side, R± being boost-invariant, the Mi-dependence of R± is independent 
of the longitudinal expansion of the source and reflects only its transverse expansion. 
Together with the breaking of the Mi-scaling of the YKP radii, the Ml- dependence of 
R± can thus be used to extract quantitatively the transverse expansion velocity of the 



source 46, 24 



4.3. Resonance decays 



As already mentioned pions in particular have the problem that in high-energy 
collisions only a fraction of about 50% or less of all pions come directly from the 
decoupling source while the rest are produced after decoupling by the decay of unstable 
resonances. The above considerations presuppose that the resonance decays do not 
affect the Mi-dependence of R±. This is not obvious, in particular since it is known that 
the Mi-dependence of the single-pion spectrum is very strongly affected by resonance 
decays [ p| . Since resonance decays contribute more to pions than to kaons they 
may also affect the M^- scaling arguments. The work by the Marburg group [ |40| on 
resonance decay effects on HBT in the context of hydro dynamical simulations indicates, 
within the standard Cartesian framework and without accounting for the cross-term, 
a possible additional Mj_-dependence of the transverse radius. However, a systematic 
analysis of resonance contributions to HBT as a function of various characteristic source 



parameters is only now becoming available [ |5(J, |5T 
4-3.1. Formalism 

Not much is known directly from experiment about the amount of resonance pro- 
duction in heavy-ion collisions. For pp collisions at similar energies a thermal model, 
where resonances are produced with thermal abundances in mutual chemical equilib- 



rium, appears to work surprisingly well [ |52], [49], [531] . We therefore assume the same to 



hold for heavy-ion collisions. We also assume for simplicity that all hadrons decouple 
at the same point. Thus all resonances are assumed to have the same emission function 
( |B0| - |82| ), adjusted only for the particle rest mass its spin J i; and its chemical poten- 
tial /ij. In chemical equilibrium /ij is given in terms of the two independent chemical 
potentials fib and fi s which account for conservation of baryon number and strangeness: 

fii = k fi b + Si fi s , (103) 

where bi,Si are the baryon number and strangeness for resonance species i. fi s is 
determined in terms of T and fib by the condition of strangeness neutrality of the 



collision region [ |54| which cannot be violated on the time scale of strong interactions. 
For illustration we will below consider the case fi^ — fi s — 0. 

I will now describe very shortly the formal steps for calculating the resonance 
contributions to the correlation function; for more details I refer the reader to Ref. [ 
50 1 . The total source distribution of pions (and similarly for other stable particle 
species, although for pions the resonance contributions are most important) can be 
written as 

S^x,p) = Sf r (x,p) + S r ^(x,p) , (104) 

where the first term on the r.h.s. is the contribution from the directly emitted pions 
(Eq. (|80| ) with J = 0, m = m n ) and the sum contains all contributions from resonance 
feed- down: 

/s + r d^P 

dsg(s) J —5(p-P-M r E*) 

J d A X J™ dr T e" rr 5 (4) (x - (x - -^-r) ) S^{X, P) . (105) 



x 



Here capital letters indicate coordinates associated with the parent resonance r, lower 
case letters are associated with the decay pion. s is the the invariant mass of the other, 
unobserved decay products; in an n-body decay, it can vary between s_ = (J2i=2 m «) 2 



and s + = (M r — m,,-) 2 . g(s) is the decay phase space for the (n — 1) unobserved particles; 
for the isotropic 2-body decay of an unpolarized resonance it is given by 

g(s) = b ^6(s-mi), (106) 

(where b r ^ n is the branching ratio for the decay channel), and for isotropic 3-body 
decays by [ [5|| 

M r b r ^ yj[s - (m 2 + m 3 ) 2 }[s - (m 2 - m 3 ) 2 ] 

2-7TS Q(M r ,m 7T ,m 2 ,m 3 ) 

f s + ds 1 i / 

Q(M r ,m n ,m2,m 3 ) = / — \J{M r + m^) 2 - s'Js + - s' 

J S— S 

x \js- — s'\J (m 2 - m 3 ) 2 - s' . 
p*, E* are the momentum and energy of the decay pion in the resonance rest frame, 

/ -— - , J[(M r + m n ) 2 - 8 ][(M r - jjg - s] 
E =sjml+p*\ p = — , (108) 



and functions of s only. The r-integration in (|105|) extends over the exponential decay 
probability of the resonance with total decay width T. The 4-dimensional ^-function 
of the space-time coordinate X ensures that for the pion to appear at point x from a 
resonance decaying at time r, the parent resonance with momentum P must have been 
emitted from the source at point X — (P/M r )r. 

The integration over the resonance momentum P is restricted by the energy- 
momentum constraint 5(p-P — M r E*). In the coordinate system where the momentum 
p of the decay pion is given by 

p^ = (m_i_ coshy,p±, 0, m± sinh y) , (109) 

(see (|8~9"D), the resonance momentum P is parametrized as 

= (Mi cosh Y, P ± cos $, P± sin $, M± sinh Y) . (110) 

For p± ^ the 5-function can be used to fix the azimuthal angle $ of the resonance 
momentum P to 

EE p -p T P T -E*M m ± M ± cosh(y - y) - E*M 

$± = ±$ with cos $ = £ — = v - * ; . 

P±P± P±P± 

(111) 

Let us denote by P the two values of P obtained by inserting the two solutions Ql 11 ) 
into ( 110 ). After doing the $- and X-integrations in ( |105| ) one thus obtains 

S r ^Jx,p) = M r / dsg(s) dY ' dM\ \ drT<r Vr 
Js_ JY- Jm 2 Jo 



lp\{M\ - M 2 ) - [E*M r - m ± M ± cosh(F - y)} 2 ' 



with the kinematic limits 

M ±± = M ± ±AM ± (113) 



E*M r m ± cosh(r - y) M r p± \J E* 2 +p\-m\ cosh 2 (F - y) 
m\ cosh 2 (y — y) — p\ w? L cosh 2 (F — y) — p\ 



Y± = y ± AY = y ± In I — + 



2 

ml 



resulting from the zeroes of the square root in the denominator (which, incidentally, 
can also be written as p±P±\ sin <3>|). - For the limiting case p± = 0, the constraint 
p-P = M r E* cannot be used to do the ^-integration. One then uses it to do the 
Mj_-integral: 



S r ^(x;y,p ± = 0) = M r f + ds g{s) f d<t> f + dY 

Js_ JO JY- 



Y + „, M r E* 



Y- w? cosh 2 (y — y) 
P 



(115) 



M 



M r E* 
ttity cosh(V — y) 



For the more generic case p± ^ a few further manipulations are useful in practice 
50f : Rewriting the square root in ( |112|) as 



1 1 



'm\ cosh 2 (F -y)-p\ \/(AM ± ) 2 - (M ± - M ± ) 2 
and introducing new integration variables v G [—1, 1], £ G [— n, tt] via 



(116) 



M± = M ± + AM ± cosC, (117) 
Y = y + vAY, (118) 

Eq. ( |112j ) can be further transformed into 

S R ^(x,p) = £ I^J o °° drTe-^Sf (x - ^rT, pA , (119) 
with the following shorthand for the integration over the resonance momenta: 

= M r dsg(s) / d([M ± + AM ± cosC) . (120) 

JR J a- J-l J m 2 rn!i h 2 (>»AY\ - r?. JO V ' 



1 y^J-Cosh (vAY) — p\ 

For the calculation of the correlation function we need the Fourier transform of the 
emission function 

S r ^(q,p) = f d 4 xe^S r ^P) =E / : \. P ± S^faP*), (121) 
and must evaluate 

c{ (g,K)| 2 + 2E^Re[^ ir (g,K)g r .^(g,K)] + |E^^(g,^)| 2 

(122) 

where the denominator includes all contributions. Numerically, this is a rather involved 
expression. If the resonance contributions are small compared to the direct term one 
can use the Grassberger approximation [ [5(| in which the last term in the numerator is 
neglected. For heavy-ion collisions this is not good enough since about 50% of all pions 
come from resonance decays. Instead, we can try to exploit the connection from Lecture 
2 between the half- widths of the correlation function and the space-time variances which 
are now given by 

/~ ~ \(TS-\ _ Z)r J d X X^Xy Sr^njx, K) (ioo\ 



Here the sum now runs over all contributions, including the direct pions. It is instructive 
to rewrite the average over the emission function in the following form: 



(x v ){K) = J2fr(K) (Xu)r(K), 
r 

{x^c v ){K) = J2M K ) (x»x v ) r {K). (124) 



Here we introduced the single-particle fractions [ |4 

f(K) Sd*xS r ^(x,K) dWjd^K Vff __ 1 , 

M } ~ ErI^xS r ^,K)- dN** ' 2-.tr{K)-l, (125) 

which give the fraction of single pions with momentum K resulting from decay channel 
r, and the average (. . .) r with the effective pion emission function arising from this 
particular channel: 

i---u*= %££io - (126) 

The variances (|123|) can then be rewritten as 

r rr i 

The first term has an easy intuitive interpretation: each resonance decay channel r 
contributes an effective emission function S r -> n . The full variance is calculated by 
weighting the variance (homogeneity length) of the emission function from a particular 
decay channel with the fraction f r with which this channel contributes to the single 
particle spectrum. The second term in (|127|) is due to the fact that in general the 
effective emission functions from the various decay channels have different saddle points; 
it somewhat spoils the intuitive interpretation of (123). 

4-3.2. Influence on HBT radii and non-Gaussian features 

It turns out that, contrary to the situations discussed before, in the case of long- 
lived resonances the expressions ( 123|) are not very useful for a quantitative under- 
standing of the correlator, although certain qualitative features can still be extracted 
relatively easily. The reason for this is best explained by considering the simple exam- 
ple of only one longlived resonance in a 1-dimensional space. Let us model the emission 
function for the direct pions by a Gaussian in coordinate space with width R^ ir and 
(somewhat unrealistically) the emission function of the pions from the decaying reso- 
nance by a second Gaussian with much larger radius -Rhaio (assuming that the resonance 
travels on average a distance of order -Rhaio before it decays), with weights e and (1 — e), 
respectively: 

S n {x, K) = Sf T (x, K) + S r ^(x, K) = (1 - e) e - x2/(2i ^) + e e-^ 2 ^) . (i 2 8) 

Then the correlator is given by 

C(q, K) - 1 = (1 - ef e~ R2 ^ 2 + £ 2 e'^ 2 + 2 e(l - e) e -^L+*LJ« 2 /i . (129) 

If e is small, but -Rhaio is large, then the correlator is a superposition of a large, broad 
Gaussian with width l/i?dir an d weight (1 — e) 2 , a second, narrower Gaussian with 
width J2/(R\ r + -Rh a i ) and smaller weight 2e(l — e), and a third, extremely narrow 



Gaussian with width l/i?haio an d tiny weight e 2 . Obviously, the rough structure of the 
correlator will be determined by the large and broad direct contribution; the two other 
contributions will, however, modify its functional form: 

(i) If the resonance is shortlived such that -Rhaio > -Rdir its effect on the correlator will 
be minor; its shape will remain roughly Gaussian, with a width somewhere between 
l/i?dir and l/i?haio, depending on the weight e of the resonance contribution. 

(ii) If the resonance lifetime and thus -Rhaio are extremely large, the second and third 
term in (|129|) will be very narrow and, due to the finite two-track resolution of every 
experiment, may escape detection; then the correlator looks again Gaussian with a 
width 1/R dil , but at q = it will not approach the value 2, but 1 + (1 — e) 2 < 2. The 
correlation appears to be incomplete, with an "chaoticity parameter" A = (1 — f r ) 2 = 
(l-e) 2 . 

(iii) If the resonance lifetime is in between such that -Rhaio -Rdir but 1/irWo being 
still large enough to be experimentally resolved, all three Gaussians contribute, and 
the full correlator deviates strongly from a single Gaussian. 

In cases (ii) and (iii) the space-time variances give misleading or outright wrong 
results for the width of the correlation function. As noted in connection with Eq. (|53|), 
they reproduce the curvature of the correlator at q = which for our toy model is 



1 d 2 C(q) 



2 8q 2 



= (l-e)R 2 div + eRl lo . (130) 

9=0 



In case (ii), for not too small values of e, this is dominated by the second term although 
the resonance contribution is not even visible in the correlator! On a quantitative level, 
the situation is not very much better for case (iii). 

For the case (ii) of very long-lived resonances there is, of course, an easy way to 
save the usefulness of the space-time variances: if one simply leaves them out from the 
sum over resonances in (|127 ), but only includes them via an "chaoticity parameter" 



X(K) = 1 - £ f r (K) , (131) 

V r=longlived / 

the roughly Gaussian contributions to the correlator from the direct pions and short- 
lived resonances are still correctly reproduced. The real head ache comes from reso- 
nances with an intermediate lifetime which lead to a large halo but can still be experi- 
mentally resolved. They cause appreciable deviations from a Gaussian behaviour for 
the correlator and cannot be reliably treated by the method of space-time variances. 

In nature there is only one such resonance: the u meson, with a lifetime of ap- 
proximately 20 fm/c. All other resonances either live so shortly (typically 1 fm/c) that 
they hardly modify the correlator, or so long that their contribution to the correlator 
cannot be resolved such that they only affect A. At low K±, however, up to 10% of all 
pions come from cu-decays (f w (K = 0) ~ 0.1), and their non-Gaussian effects on the 
correlator can be clearly seen. On account of the u a full numerical evaluation of the 
correlator ( |122j ) or a treatment with more powerful analytical methods which can deal 
with the non-Gaussian features of the correlator (q- variances, see Ref. [ [50|D become 
indispensible. 

4-3.3. K -dependence of correlator including resonance decays 



I now show some numerical results for the correlation functions resulting from the 
emission function (|80|), but now including the resonance contributions. The complete 



spectrum of relevant resonances is included, and in the decays the 2- and 3-body decay 
kinematics is fully taken into account. The HBT radii are extracted from a Gaussian 
fit to the numerically calculated correlation function. A detailed technical discussion is 
given in Ref. [ (50"|. 
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Fig. 7. The influence of resonance decays on the A/^-dependence of R s (a,b) and R a 
(c,d) for Y cm = pion pairs. a,c: no transverse flow; b,d: transverse flow rapidity 
i]f = 0.3. The Gaussian transverse radius is here R — 5 fm, and T = 150 MeV. The 
HBT radii are extracted from unidirectional fits to the correlator in the respective 
direction of q. (Figure taken from Ref. [ |50|.) 

Fig. 7 shows results for the standard Cartesian parameters R s and R Q from 1- 
dimensional fits to the numerically computed correlator in the respective (q s or q ) 
directions (setting the other components of q to zero). One sees that the effects of the 
short-lived resonances with lifetimes of order 1 fm/c on R s are essentially negligible, 
both at vanishing and at nonzero transverse flow. Only the uj with its intermediate 
lifetime of 20 fm/c affects R s , but only for vanishing transverse flow. There it induces a 
weak M^-dependence at small M±_ even in the absence of transverse flow; at Mj_ > 500 
MeV the contribution of the u dies out, and R s again becomes M^-independent (which 
would not be the case if it were affected by flow). At rjf = 0.3 and 0.6 [ |5(J not 
even the uj generates any additional Mi-dependence! -R D shows some effects from the 
additional lifetime of the resonances, in particular from the long-lived uj. Resonances 
with much longer lifetimes than the u (in particular all weak decays) cannot be resolved 



experimentally in the correlator and have no effect on the HBT radii. 
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Fig. 8. Same as Fig. 7, but now with parameters extracted from a complete 5- 
dimcnsional fit to the correlator (see text). Top row: R s ; second row: R Q ; third row: 
Ri; bottom row: effective chaoticity parameter A. Left column: no transverse flow; 
right column: transverse expansion with rjf = 0.3. The cross term R Q i vanishes for 
Y = pairs. (Figure taken from Ref. [ |5Q].) 

The weak effect of resonances on R s = R± seems surprising: due to their non- 
zero lifetime they should be able to propagate outside the original source before decay 
and form a pion "halo" [ [57|. This effect is, however, much weaker than naively 
expected: most of the resonances are not very fast, and the halo thickness is thus only 
a fraction of the resonance lifetime. At finite transverse flow an additional effect comes 
into play: it turns out that then the effective size of the emission function for directly 
emitted resonances is smaller than that for direct pions [ [5(J]! At 7^=0.3 and 0.6 this 
even slightly overcompensates the halo effect, and altogether the resonances change 
neither the size nor the Mi-dependence of R s . 

In Fig. 8 also the longitudinal radius Ri is shown. Here even the shortlived reso- 
nances are seen to make an effect. It can be essentially traced to a lifetime effect: Since 
the pions from the short-lived resonances appear typically 1 fm/c later than the direct 



pions, the source has in the meantime expanded longitudinally to a situation with a 
smaller longitudinal velocity gradient (the latter goes like 1/r). Thus the resulting pion 
source has a larger longitudinal length of homogeneity and features a larger value Ri. 



4-3.4- The "chaoticity parameter" A 

As noted above, the c^-decays make the correlator non-Gaussian. Slight non- 
Gaussian features exist even without resonance decays in the longitudinal direction, 
induced by the strong longitudinal flow, and for large transverse flow rjf become also 
visible in the transverse direction. A Gaussian fit to such a slightly non-Gaussian cor- 
relator in general does not extrapolate to the correct value at q = 0, but introduces 
an effective "chaoticity parameter" A. In 1-dimensional fits along the three cartesian 
directions of q the different degrees of non-Gaussicity lead to different values of A. This 
situation is aggravated by the resonance contributions, in particular from the u which 
affects the correlator differently in each direction. Still a different value of A is found 
in a 5-dimensional fit to the correlator, using the 4 standard Cartesian HBT radii and 
A as fit parameters, because now A has to "compromise" between the values found in 
the unidirectional fits. From the last row in Figs. 8 one sees that even without the 
contribution from the uj and the longlived resonances the effective A is below 1 by up 
to 10%; inclusion of the uj reduces it further to about 80-85% at K = 0. Thus not only 
the very longlived resonances affect A as anticipated above, but so do to some extend 
the medium- and shortlived resonances and even flow. 

Comparing Figs. 7 and Fig. 8 one sees that this compromise in A between the 3 
different unidirectional fit values and the one resulting in the 5-dimensional fit also 
changes the fitted HBT radii. On the same (small) level as already observed for the 
contribution from uj decays, it even affects the M^-dependence of R±. Both effects, 
the non-Gaussian features of the correlator and the effective "chaoticity parameter" 
which differs from A = 1, vanish at large K because the resonance decay pions are 



concentrated at low K [ 49] 



5. FINAL REMARKS 

In these three lectures I have presented to you two-particle intensity interferometry 
for relativistic heavy-ion collisions both as an intellectually stimulating problem and 
as a powerful practical method. Its application to nuclear collisions has turned out 
to be much more difficult than expected from the astrophysical analogue (two-photon 
intensity interferometry of stars). But at the same time, due to the dynamical nature 
of the problem, the physics of heavy-ion collisions is very much richer, and I have tried 
to show you that the HBT method is up to the task of clarifying a lot of this physics 
in a rather direct manner. 

The key to our understanding of the HBT method and how to apply it to dynamical 
situations are the model-independent expressions derived in the second lecture, which 
express the HBT width parameters in terms of second order space-time variances of 
the emission function. They provide the basis of a detailed physical interpretation of 
the measured HBT radii. They show that generally the HBT radius parameters do 
not measure the full geometric extension of the source, but regions of homogeneity 
inside the effective emission function for particles with certain fixed momenta. For 
expanding systems these are usually smaller than the naive geometric source size and 
decreasing functions of the pair momentum. For systems with finite lifetime the HBT 
parameters usually mix the spatial and temporal structure of the source, and their 



unfolding requires model studies. 

With the new YKP parametrization a method has been found which, for systems 
with dominant longitudinal expansion, cleanly factorises the longitudinal and transverse 
spatial from the temporal homogeneity length. The effective source lifetime is directly 
fitted by the parameter R ; it is generically a function of the pair momentum and largest 
for pairs which are slow in the CMS. Another fit parameter, the YK velocity, measures 
directly the longitudinal velocity of the emitting fluid element, and its dependence on 
the pair rapidity allows for a direct determination of the longitudinal expansion of the 
source. Without transverse expansion, the YKP radius parameters show exact M±- 
scaling. The breaking of this scaling and the M±- dependence of the transverse radius 
parameter R± allow for a determination of the transverse expansion velocity of the 
source. Resonance decays were shown to mostly affect the lifetime parameter and, as a 
consequence, the longitudinal homogeneity length. They leave the Mi-dependence of 
R± nearly unchanged and thus do not endanger the extraction of the transverse flow 
via HBT. 

With this new and detailed understanding of the method, I believe that HBT 
interferometry has a begun a new and vigorous life as a powerful tool for reconstructing 
the geometric and dynamic space-time characteristics of the collision zone from the 
measured momentum spectra. 
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Keywords: Two-particle interferometry / Hanbury-Brown-Twiss interferometry / 
HBT radii / relativistic heavy-ion collisions / quark-gluon plasma / resonance decays 
/ collective flow / hydrodynamics / YKP parametrization / Mperp-scaling / emission 
function 

Abstract: I discuss two-particle intensity interferometry as a method to extract from 
measured 1- and 2-particle momentum spectra information on the space-time geometry 
and dynamics of the particle emitting source. Particular attention is given to the rapid 
expansion and short lifetime of the sources created in relativistic heavy-ion collisions. 
Model-independent expressions for the HBT size parameters in terms of the space- 
time variances of the source are derived, and a new parametrization of the correlation 
function is suggested which allows to separate the transverse, longitudinal and tempo- 
ral extension of the source and to measure its transverse and longitudinal expansion 
velocity. The effects of resonance decays are also discussed. 



